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APPROXIMATE METHOD FOR CALCULATING TRANSONIC FLOW 
ABOUT LIFTING WING- BODY CONFIGURATIONS 

Richard W. Barnwell 
Langley Research Center 

SUMMARY 

The three-dimensional problem of transonic flow about lifting wing -body configu- 
rations is reduced to a two-variable computational problem with the method of matched 
asymptotic expansions. The computational problem is solved With the method of relaxa- 
tion. The method accounts for leading-edge separation, the presence of shock waves, 
and the presence of solid, slotted, or porous tunnel walls. The Mach number range of 
the method extends from zero to the supersonic value at which the wing leading edge 
becomes sonic. A modified form of the transonic area rule which accounts for the effect 
of lift is developed. This effect is explained from simple physical considerations. 

INTRODUCTION 

It is well known that there are approximate methods for calculating transonic flow 
about lifting swept-wing configurations with span-length ratios of order one and small 
reduced span-length ratios ^the reduced span is the product of the span and the factor 

\/|l " ^00^1 where M^ is the free -stream Mach number^ if the angle of attack is suffi- 
ciently small. A recent example of such a method is that of Stahara and Spreiter (ref. 1). 
These methods are based on the analysis of Heaslet and Spreiter (ref. 2) which shows 
that if the angle of attack is less than the wing thickness ratio, the problem of transonic 
flow over a wing can be separated into noninteracting thickness and cross-flow problems. 
The analysis also shows that at these small angles of attack, the effect of lift on the flow 
in the outer region is negligible and that the determination of the shock-wave strength 
and location is part of the thickness problem. It should be noted that these results are 
consistent with those of Oswatitsch and Keune (ref. 3), who have proved that at transonic 
speeds the flow in the outer region of a nonlifting wing is mathematically equivalent to 
that about the axisymmetric body with the same cross-sectional area distribution (the 
equivalent body), and Whitcomb (ref. 4), who has demonstrated experimentally that the 
transonic wave drag of a nonlifting wing -body is almost the same as that on the equivalent 
body. The basic solution to the cross-flow problem is the slender-wing solution of 



Jones (ref. 5). This solution has been generalized for lifting wing -body configurations 
by Spreiter (ref. 6) and Ward (ref. 7). It should be noted that slender-wing theory is 
applicable to wings with span-length ratios of order one only for Mach numbers very 
close to one and not throughout the entire transonic range. 

Panel methods based on the linearized potential equation have been applied to 
lifting configurations in the high subsonic and low supersonic flow regimes where tran- 
sonic phenomena are present. However, these methods cannot be applied throughout the 
entire transonic range because of the obvious inability of linear methods to resolve the 
shock waves associated with transonic flow properly. 

At the present time three-dimensional finite -difference methods are being developed 
for application to the problem of transonic flow about lifting wing-body configurations. 
These methods, which are formulated in terms of the velocity potential, employ the suc- 
cessive line overrelaxation procedure introduced by Murman and Cole (ref. 8). Bailey 
and Ballhaus (ref. 9) and Schmidt, Rohlfs, and Vanino (ref. 10) have obtained solutions for 
infinite -cylinder swept-wing combinations and Klunker and Newman (ref. 11) have treated 
infinite -cylinder rectangular -wing combinations. It is probable that methods of this type 
will be generalized so that configurations with finite-length bodies can be treated. It 
should be noted that these methods require large amounts of computing time and com- 
puter storage since the problem being calculated is three dimensional. Also, these 
methods have not been developed to the point where the phenomenon of leading-edge sep- 
aration can be treated. 

The purpose of this paper is to present an approximate method for calculating tran- 
sonic flow about lifting configurations with span-length ratios of order one and small 
reduced-span — length ratios which are at angles of attack large enough to insure that the 
effects of lift and thickness are comparable in the outer region. The approximations 
used in the method are based on the analyses of Barnwell (refs. 12 and 13) and Cheng and 
Hafez (refs. 14 and 15), which show that in this angle -of -attack range the thickness prob- 
lem is dependent on the lift solution. These analyses are used to reduce the three- 
dimensional problem of transonic flow about a lifting swept-wing configuration to a two- 
variable computational problem which is solved with a finite -difference method of the type 
introduced by Murman and Cole (ref. 8). It should be noted that a preliminary version 
of the present method was presented in reference 16, and that some results obtained with 
the final version are given in reference 17. A description of the computer program for 
the present method and a user^s manual are given in reference 18. 
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SYMBOLS 


Where applicable, values are given in both SI and U.S. Customary Units. 

A coefficient of quasiconical lift potential (see eq. (26)) 

■^k’^k’^k coefficients in finite -difference formula (F7) given by equations (F8) 

Aj^,Ag wing-geometry parameters 

AQjAjjAg, . . .,Aj^ coefficients of Fourier series (eq. (D19)) 

a constant of proportionality in equation (B21) for angle of attack 

a ,b slot width and slot spacing 

w w 

a ,c slot-width parameter and slot- thickness parameter 

w w 

B coefficient of proportionality in equation (74) for effective cross-sectional 

area of lifting configuration 

coefficient of wall-induced contribution to velocity potential (see eq. (34)) 
Bj^jBg wing-geometry parameters 

b wing semispan 

bj^jbg apparent semispans in quasiconical model (see fig. 39) 

C body geometry parameter 

coefficient 

Cp pressure coefficient 

coefficient of wall -induced contribution to velocity potential (see eq. (35)) 
c(y) chord at span station y 
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Cj^,C 2 distances depicted in figure 39 

E(k) complete elliptic integral of second kind with modulus k 

F(x,r,o)) forcing function in equation (F2) 


Fg(x) 


f(x) 


nondimensional radius of equivalent body 
nondimensional dipole strength distribution 


fy(x) contribution to f(x) due to vortex (see eq. (27)) 

G(x,r,a)) function in equation (F3) 

Gj^(x),Gg(x) arbitrary fimctions in equations (B26) and (B28) 
g(x) lift function defined by equation (D13) 

gg(x),g 2 functions of x in inner potentials <^2 2’ ^2 1 

H(x) source -strength function defined by equation (B17) 

H(a,r) function defined by equation (D18) 

K constant of proportionality relating and 6 (see eqs. (B5) and (B6)) 

K(k) complete elliptic integral of first kind with modulus k 

k wing-geometry parameter defined by equation (C12) 

kj^,k 2 wing-geometry parameters 

k(x) lift function defined by equation (D12) 

L(x) total lift on portion of configuration between nose and longitudinal location x 

given by equation (64) 

Lq(x) lift on portion of configuration between nose and longitudinal location x due 

to nonlinear effects 
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L^(x) 


S. 

M 




m(x,y) 

mp(x,y/Y2) 

m^(x,y) 

n 


P 


R 


Rn 


Ri 


r 


rb(x) 

^0 

S(x) 


lift on portion of configuration between nose and longitudinal location x 
due to linear effects 

longitudinal length scale, configuration length unless otherwise noted 
local Mach number 

Mach number of component of upstream flow normal to shock 
free -stream Mach number 
= |(m^ - m.) 

nondimensional function in equation (B31) 
function in equation (B12) for potential near configuration 
exponent in equations for equivalent body radius (see appendix A) 
coefficient in equation (33) for porous-wall boundary condition 
function defined by equation (14) or (El) 
coefficient of equation (D29) given by equation (D30) 

coefficient of equation (D40) given by equation (D41) unless otherwise noted 

radial coordinate 

outer radial variable ur 

radius of equivalent body 

value of r at which body boundary condition is applied 
Mirels^ S function for wings with swept trailing edges 
effective cross-sectional area of lifting wing body (see eq. (74)) 
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*ref 


U„ 



V. 


W(x,X) 




w 


X 

Xy 

X 

^0 

Xi,X2,Xg 

xpXn,Xin 

Y(x,X) 

Yv(x) 


reference area for lift coefficient 
maximum radius of equivalent body 
free-stream velocity 

perturbation velocity components in axial, radial, and cross-flow directions, 
respectively 

x-component of perturbation velocity on surface of wing alone 
average velocity at points near axis given by equations (57) 

9y dz 
complex potential 

complex potential for flow due to vortices 
downwash for wing alone 

complex variable in cross-flow plane, y + iz 
value of X where vortex core is located 
longitudinal coordinate 

quantity used in expression for equivalent body (see appendix A) 
value of X where wing leading edge intersects body surface 
values of x used in description of wing planform (see appendix A) 
values of x where different segments of wing planform intersect 
complex function defined by equation (C4) 
value of Y at vortex-core location 
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y 


spanwise variable 


y^(x) spanwise location of vortex core in physical plane 

y 2 (x),y^(x) functions giving location of wing leading and trailing edges, respectively 
y y(x) , y j j(x) , y straight-line segments used to construct planform (see appendix A) 


yio’^no 


values of y where different segments of wing planform intersect 


Z complex function defined by equation (C9) 


z vertical coordinate 


value of z where vortex core is located 
a angle of attack 



j3(a) function defined by equation (D16) 

y ratio of specific heats 

vortex core strength 

A parameter used to specify wing planform shape (see appendix A) 

AF forcing function due to wall effects given by equation (44) 

Ax^,Ax 2 ,Ax 3 mesh spacings given by equations (F4) 

A 0 ^’ difference in the values 0 ^^ above and below wing 

6 equivalent -body thickness ratio, t/i 

61 , 62,63 gage functions of outer expansion (eq. (B3)) given by equations (B5) and (B22) 

77 spanwise coordinate in Z-plane given by equation (25) 
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A 

X 

>-v(x) 


a,T 


7V(x) 






values of 77 for leading and trailing edges, respectively, given by 
equations (Cll) 

polar angle of body -oriented coordinate system 
coefficient in equation (29) for slotted-wall boundary condition 
coefficient of derivative in equation (42) given by eqviation (43) 

quantity defined by equation (49) or (50) 

spanwise location of vortex core in Y-plane, see equation (C20) 

radial stretching parameter in outer region 

vertical coordinate in Z -plane 

angular coordinates given by equations (D15) 

vertical location of vortex core in Y-plane, see equation (C20) 

part of due to vortex, see equation (D33) 

second-order part of potential cp in inner region, see equation (45) 


potentials in outer expansion (eq. (B3)) for cp 


<P 


perturbation velocity potential 


coefficients of Fourier components of cp (see eq. (10)) 


0 

K 

4>a 


potential defined by equation (Dll) 
attached-flow portion of (see eq. (D33)) 

linear lift potential 


2’ ^2 l’^2 potentials in inner expansion (eq. (B7)) for <p 
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solutions to the two-dimensional Laplace equation in cross-flow plane 
CO polar angle of wind-oriented coordinate system 

V2 two-dimensional Laplace operator in cross-flow plane (see eq. (4)) 

Superscripts: 

C complementary solution 

P particular solution 

* complex conjugate 

Subscripts: 

b backward difference expression 

c centered difference expression 

j index for x-coordinate 

k index for r-coordinate 

H lower surface of wing 

le leading edge of wing 

max maximum value 

min minimum value 

t trailing edge 

u upper surface of wing 


Primes denote differentiation with respect to x. 


PROBLEM DESCRIPTION 


The configurations which are considered in this report are characterized by the 
parameters 


5.0(1) 

6 = t « 1 > (1) 

a 

Of « 1 


where b, !L, and t are the wing semispan, longitudinal length scale (configuration 
length unless otherwise noted), and maximum equivalent- body radius, and 6 and o; are 
the thickness ratio and the angle of attack in radians. The free -stream Mach number M^^ 
is restricted to the range from zero to low supersonic. In this report it is demonstrated 
that flow can be calculated about lifting configurations in free air and in solid-, slotted-, 
and porous-wall wind tunnels. 

Wind-oriented cylindrical polar and Cartesian coordinate systems are used in this 
report. These coordinate systems are shown in figure 1. Wind orientation is chosen in 
order to simplify the governing equation and to facilitate the enforcement of boundary 
conditions at the wind-tunnel walls. The origin is at .the configuration nose, the positive 
x-axis is in the free-stream direction, the positive y-axis is to starboard in the hori- 
zontal plane, and the positive z-axis is in the upward vertical direction. 


r 



Figure 1.- Wind-oriented coordinate systems. 

In this report, wing-body combinations are approximated by flat -plate wings and 
equivalent axisymmetric bodies. It is assumed that the wing plane passes through the 
axis of the body so that the angles of attack of the body and wing are the same. The 
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geometry of an approximate wing -body combination is shown in figure 2. The equivalent- 
body surface is specified by the equation 

v^{x) = tFg(x) (2) 

where Fg(x) is the nondimensional thickness distribution of the equivalent body. The 
leading and trailing edges of the wing planform are designated by the functions y2(x) 
and yj^(x), respectively. The functions F^(x), y 2(^)7 ^1^^^ which are used in this 

report are discussed in appendix A. 



It is assumed that the flow fields under consideration are isentropic. This assump- 
tion is not strictly valid when shock waves are present in the flow field. However, it can 
be shown (ref. 19) that the largest term which is affected by the nonisentropic condition 
is of the order 


(1 - V)' 


where is the Mach number of the component of the 

upstream flow normal to the shock. Terms of this order are negligible in the present 
approximation; therefore, the isentropic assumption is valid. 


Since the flow is assumed to be isentropic, the problem is formulated in terms of 
the perturbation velocity potential cp. The governing equation for cp is 


(1 - ^ 


1) 

9x 0 x 2 


+ 2 


a ip 


a (O 


2 diP 


a (p 


dr dx ar 


r>2 aoi ax acj 


(r-i)^V2 


+ / ^ 9^ (P ^ _l_ / ac> \^ a^ ^ 2 9 c a ^ (p 2 dip diP a ^ <.0 2 dip dip a ^ ip 

[dr J 0^2 p4\acu/ q^ 2 dx 3r ax dr y2 dx doj dx du) 9 r aco 
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where is the Laplace operator in planes normal to the x-axis and is written as 


v^ = — 1 9 I 1 

2 gj.2 r 9r i*2 


( 4 ) 


The approximate boundary conditions for the configuration are written as 


d(p 

8r 


drv 


^x,r|g,co^ = -sin a sin co + cos a 


(0 s o) £ 27 t) (5) 


and 


i |^(x,r,0) = -i |^(x,r,7T) = -sin a (y 2 W = ^ = yi^ 

In this report the potential cp is made nondime ns ional with where is the 

free -stream speed, and all lengths are made nondimensional with £. 

DESCRIPTION OF METHOD 


At present, nonlinear solutions to transonic flow problems can be obtained only 
with numerical methods. Numerical computation of three-dimensional mixed-flow prob- 
lems is a costly task with current methods. In this section the analytical treatment of 
reference 13, which is reviewed in appendix B, is used to reduce the three-dimensional 
problem of transonic flow about lifting configurations with aspect ratios of order one to 
two simpler computational problems. It is shown that the dominant part of the lift 
potential is governed by a linear equation; and that once this part of the lift potential is 
known, the dependence of the perturbation velocity potential (p on the cross-flow 
angle w can be determined. Consequently, the perturbation velocity potential can be 
obtained as the solution of a two-variable problem in a plane of constant w. 


Approximate Governing Equations 

From the equations in appendix B which govern the outer potentials 
and and the first- and second-order inner potentials, it can be seen that to second 
order in sin a and 6, the equation which governs the perturbation velocity potential <p 
in both the inner and outer regions at near-sonic speeds is 


1 - - (y + 1)M, 


2 Qcp 
^ dX 


^ + v^cp = 2M^^ sin^a 
3x2 


+J_ 

9r ax dr r2 9co ax aco. 




(7) 
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where 4>^ is the linear lift potential. At near-sonic speeds 6^ is the same as 
and, as a result, is governed by the equation 


= 0 (8) 

2 

The nonlinear terms in equation (7) have been made proportional to to be consis- 

tent with other transonic theories and to insure the proper behavior in the incompressible 
limit. The terms on the right-hand side of equation (7) are of higher order than necessary 
in the outer region but are included to facilitate transition between the two regions. 

For Mach numbers appreciably different from 1, equation (7) is more than adequate 
for describing the behavior of cp. However, at these Mach numbers the linear lift 
potential 0^ is governed by the equation 




3^0 


3x' 


r * = ° 


( 9 ) 


rather than by equation (8). 


Solutions for Linear Lift Potential 

At near-sonic speeds the linear lift potential 0^ is governed by equation (8) and, 
as a result, is obtained from slender-wing theory. The Brown and Michael model 
(ref. 20) for leading-edge separation is used at these speeds. Slender-wing solutions 
for 0^ for attached and separated flow are given in appendix C. 

The linear lift potential 0^ is governed by equation (9) at subsonic and supersonic 
speeds. In this report 0^ is obtained from the theory of Lawrence and Flax (ref. 21) 
for subsonic flow and from a modified form of the theory of Carafoli et al. (ref. 22) for 
supersonic flow. A modified form of the Brown and Michael model for leading-edge sep- 
aration which accounts for Mach number effects is used for subsonic and supersonic 
flows. These methods for determining 0^ for attached and separated flows at subsonic 
and supersonic speeds are presented in appendix D. 

The theory of Lawrence and Flax for subsonic flow and the modified theory of 
Carafoli et al. for supersonic flow were chosen for use in this report because they both 
reduce to slender-wing theory in the limit as the free-stream Mach number M^ 
approaches 1, and because they are simple theories which do not involve much numerical 
computation. Consequently, the linear lift solutions which are used in this report vary 
consistently with Mach number in the range from incompressible to low supersonic flow 
and are easy to implement. However, it should be noted that more sophisticated subsonic 
and supersonic linear lift theories could have been employed for both attached and sepa- 
rated flow. 
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Two-Variable Computation of Perturbation Velocity Potential 

The perturbation velocity potential cp is calculated with an approximate equation 
obtained by evaluation of the right-hand side of equation (7) and the derivative 
in that equation with the known solution of the linear lift potential (f>^. This approximate 
equation, which has the two independent variables x and r, is solved with the method of 
relaxation. Different approximations for are used in the outside flow field, where 

r > y 2 , and the inside flow field, where r < yg. 

Outside flow field.- At transonic speeds the outside flow field ^r > includes the 
outer part of the inner region (where the governing equation is parabolic), the outer region 
(where the governing equation is of mixed elliptic -hyperbolic type), and the far field (where 
the governing equation is elliptic if the free-stream Mach number is subsonic and hyper- 
bolic if the free-stream Mach number is supersonic). From equation (B20) or equa- 
tions (B26), (B27), and (B28) it is seen that the perturbation velocity potential <p and the 
linear lift potential for configurations with one plane of symmetry can be expressed 

in the forms 


p 

(p{x,r,u)) = (Pq{x,t) + sin a ^j(x,r) sin oj + sin^of 2co + 


( 10 ) 


and 


0Q, = <^i(x,r) sin oj + . . . 


( 11 ) 


where <Pq is of second order in sin a and 6. It is tempting to separate equation (7) 
into its Fourier components and solve the component equations simultaneously. However, 
this procedure cannot be used for transonic flow fields because the basic nonlinear term 


1 -mJ -(y + 


a>", which depends on the angle 


o), would be represented by the 


term 


mJ - (y + 1)M^2 ^o’ 


I^q", which is independent of w. As a result, the 

Fourier analysis would indicate that the locations of shock waves are independent of u>, 
whereas it is clear that, in general, the shock -wave locations depend on o). 


Another approach, which is the one employed in this paper, is to use equations (10) 
and (11) to approximate the derivative and the terms involving <p^ in equa- 

tion (7) and to integrate the resulting equation in half planes of constant co. With this 
method the shock-wave locations can vary with w, and the computational problem is 
reduced to a two-variable calculation in the half planes. A separate integration must be 
performed for each value of w for which results are desired. In the horizontal plane 
(co = 0), results are calculated both above and beneath the wing. 
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The functions (jcij^(x,r) and cp 2 {x,r) which are used to evaluate the derivative 
d^(p/du)^ are approximated analytically in this paper. Account is taken of the fact that 
the forms of these functions are different for large and small values of r since the 
governing equation is parabolic near the configuration and either elliptic or hyperbolic in 
the far field. For the present it will be assumed that the lifting configuration is traveling 
in free air. 

Consider the function cp-^ for <1. It can be shown from the results of appen- 
dix C that in the parabolic inner region where slender-wing theory applies, the form is 

<^l(x,r) = M (12) 


Klunker (ref. 23) shows that the far- field solution for is 

<?l(x,r) = A(l+|)f(l) (13) 

where 



Note that the x-coordinate is measured from the nose of the configuration, and that the 
configuration has a nondimensional length of 1. In this paper the function cp^ is approxi- 
mated in the outer part of the flow field by the composite solution 


a>l(x,r) = ^ 


12 I f (x) 



(15) 


for M <1, where 

f(x) = f(l) (x > 1) 

f(x) = 0 (x < 0) 


Note that this equation is correct in the asymptotic limits of the parabolic inner region, 
where equation (12) applies, and the elliptic far field, where equation (13) applies. It is 
assumed to be correct in the intervening mixed outer region. 
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Consider the function <p^ for >1. In appendix E it is shown that in the far 
field in the interval 


0 ^ X - /Sr £ 1 


( 16 ) 


the function <p^ can be approximated as 




r(l - /Sr/x)2 


\/l - ^ 




(17) 


where 




/3=\/M^ -1 


(18) 


In this report the function cp.^ is evaluated with equation (17) in the outer part of the 
flow field in the interval given by inequalities (16). Outside of this interval the func- 
tion cp.^ is assumed to vanish. In order to obtain equation (17), it is necessary to 
assume that the variation of y 2 with x is essentially linear. It should be noted that 
the right side of equation (17) approaches that of equation (12) as r becomes small. 

Consider the function (^ 2 * It can be seen from equation (B20) that this function 
can be written in the inner region as 




(19) 


Note that the first term in equation (19) does not approach zero as r becomes large. 
In order to insure the proper asymptotic behavior for large values of r, it is assumed 
that for < 1, the function <p 2 is given by the equation 




( 20 ) 


in the outer part of the flow field. For > 1, cp 2 is evaluated with equation (19) in 
the region 


r < x//3 
0 £ X S 1 

Outside of this region </?2 is assumed to vanish. 


( 21 ) 
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It is shown in appendix C that the function f(x) for attached transonic flow past 
configurations with straight trailing edges is 



1 + 




(22) 


In calculating the derivatives f and f", it is assumed that r^^* is negligible compared 
with Yg'. Although the quantity (r] 3 /y 2 )^ is very small except near the point where 
the wing leading edge intersects the body, it is included so that the derivative f'(x) will 
vanish at the point where the wing leading edge intersects the body rather than at the 
apparent wing apex inside the body. The derivative f*(x) for values of x where the 
wing has both a leading edge Y 2 (x) and a trailing edge yj(x) can be evaluated exactly 
with slender-wing theory. However, this evaluation is very tedious. In this report an 
approximate form, which is a generalization of the approximate form developed for wings 
by Mangier (ref. 24), is used. This approximate form is written as 


|=0.85y2 




y2 -yi 

F2 + yi 


(23) 


The function f(x) is determined by numerical integration of equation (23). 

For subsonic Mach numbers appreciably different from 1, solutions to equation (9) 
for for attached flow are obtained with the approximate method of Lawrence and 

Flax (ref. 21) for wing-body configurations. A brief discussion of this method, which is 
an extension of that of Lawrence (ref. 25) for wings alone, is given in appendix D. The 
method is applicable to configurations with unswept trailing edges and aspect ratios of 
order one or less. With this method solutions for the function f(x) are obtained which 
to lowest order satisfy the trailing-edge Kutta condition. The terms in equations (12), 
(13), (15), (17), (19), and (20) which depend on the function f(x) and its derivatives are 
evaluated with this solution. It should be noted that the slender-wing solution for f(x) 
does not, in general, satisfy the Kutta condition at the trailing edge. 

It is shown in appendix D that the function f(x) is related to the values of cp^ on 
the upper and lower surfaces of the wing, which are designated as <p^ ^ and cp^ 
respectively, by the equation 


11 




(24) 
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where the transformed coordinate -q is related to the physical coordinates y and z 
by the equation 


q 



(25) 


Since the function f(x) is determined directly with the method of Lawrence and Flax, an 
additional assumption concerning the spanwise variation of (j>^ must be made in order 
to evaluate the pressure distribution. In this report it is assumed that the spanwise vari- 
ation of on the surface is the same as that predicted by slender -wing theory. 

For supersonic attached flow an approximate solution to equation (9) is obtained 
with a method which is based on the quasiconical method of Carafoli et al. (ref. 22) for 
supersonic flow past wings with subsonic leading edges. The supersonic applications of 
the present method and the method of reference 22 are restricted to wings with unswept 
trailing edges. It is assumed that the forms of the potential <f)^ on the upper and lower 
surfaces of the wing, respectively, can be approximated as 




(26) 


where q is given by equation (25). The quantities A, bj^, and b 2 are slowly varying 
functions of x and y (or z) which correspond to the quantities A*, and in 
reference 22. The method which is used to evaluate the functions A, bj, and b 2 is 
given in appendix D. 

The function f(x) for attached flow is obtained from equation (24) by numerical 
integration. Values for the integrand are obtained from equation (26). The derivatives 
of f(x) used in equations (19) and (20) are determined numerically with upwind difference 
formulas. 

If the flow is separated at the leading edge of the wing so that there is a vortex 
above the wing near the leading edge, or if the flow is attached at the leading edge but 
there is a vortex over the wing which was generated by a strake ahead of the wing, there 
is an additional contribution to the function f(x) of the form 


f (x) = 2y^(x) X, fx) 


Y2(x) + 


-b W 


Y2(x) 


( 27 ) 


An expression for the product yy(x)Xy(x) can be obtained from equation (C26). 
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Consider the case of flow past a model in a tunnel. If the tunnel walls are solid, 
the boundary condition at the wall (r = ^max) 


ar 





( 28 ) 


Baldwin, Turner, and Knechtel (ref. 26) show that if the tunnel has slotted walls, the 
boundary condition is 


I.l^- 0 

ax ar K ax 


= ^max) 


(29) 


where k is related to the slot width a^ and the slot spacing b^ by the equation 





(30) 


Chen and Mears (ref. 27) also obtain equation (29) as the boundary condition for slotted- 
wall tunnels where k is related to the slot-width parameter a^, slot spacing b^, 
and slot- thickness parameter c^ by the equation 



It should be noted that no assumptions concerning Mach number were made in the deriva- 
tion of equation (29). In this report it is assumed that the slot geometry does not vary 
with X. Consequently, equation (29) can be integrated with respect to x to yield 

This form is used in this report. For tvinnels with porous walls the boundary condition 
is 


ar 9x 


= ^max) 


(33) 
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It is assumed that the tunnel wall lies in the region where the Fourier expansion 
(eq. (10)) applies. Since the boundary conditions are linear, each component on the right 
side of equation (10) must satisfy them. Assume that the functions cp-^{x,r) and ^ 2 (^ 5 ^^) 
are of the form 

^j(x,r) = + B^(x)r (34) 

^2(x.r) = ('£" - + C„(x)r2 (35) 

where the coefficients B^(x) and C^(x) are to be determined from the boundary con- 
ditions. The forms of the last terms in equations (34) and (35) are chosen so that the 
terms in equation (10) proportional to B^(x) and C^(x) are solutions to the two- 
dimensional Laplace equation. 

When equations (34) and (35) are substituted into the slotted-wall boundary condi- 
tion (eq. (32)), it can be seen that the quantities B^(x) and C^(x) are given by the 
equations 


and 


B^(x) 


f(x) 

^Lax 


^ - ^maxA 

^ ^max/^ 


Cw(x) 


4f(x) f'(x)(2 - 

) + (y + l)f'(x) f"(x) r^^ 

«dax( 

2 + r /fc^ 

max/ j 

1 


(36) 


(37) 


Note that the solid-wall boundary condition (eq. (28)) is simply the special case of the 
slotted-wall condition (eq. (32)) with l//c = 0. 

When equations (34) and (35) are substituted into the porous-wall boundary condition 
(eq. (33)), it is found that B^(x) and C^(x) are governed by the ordinary differential 
equations 




max 


“ o 


f’ - 


Pr 


max 


max 


(38) 
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and 




-Cw = 

^*^max ^^max 


(ff)' 


Pr 


-ff' 
max J 


(39) 


The solutions to equations (38) and (39) are 


exp 


(- " 


Bw(x) = - 


\ ^""max/ (’S=x 


max 


) exp 
S=-oo 


/ 2s ^ 

d 

y^^max^ 

Ids 


f(s) exp /- 

I Pr 


max 


ds 


(40) 


and 



To summarize, the equation which is solved in the outside flow field is 


A = -^[sin a (/7,(x,r) sin o) + 4 sin^o! (p^(x,r) cos 2u) 

dx^ r 8r \ 8r / r^L 


^ 2M 2 sin^a 

oo A 


where 


A = 1 - - (y + mj 

oX 


(42) 


(43) 


and where cp-^ and <p 2 are given by equations (15) and (20) for g 1 and equa- 
tions (17) and (19) for >1. K a timnel wall is present, equation (28), (32), or (33) 
is used for the boundary condition, and the quantity 
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AF = i sin a B^(x) sin o) + 4 sin^O! C^(x) cos 2 cd 


(44) 


is added to the right side of equation (42). The functions B^(x) and C^(x) are 
obtained from equations (36) and (37) for solid or slotted walls and equations (40) and (41) 
for porous walls. 

Flow field near configuration.- Consider the flow in the close vicinity of the config- 
uration. Assume that in this region the perturbation velocity potential cp can be 
approximated as 


^(x,r,w) = sin a + 4>Q(x,r,w) 


(45) 


where the potential <p^ is the linear lift potential and the potential «&q is of second 
order in the parameters 6 and sin o!. Solutions for cf>^ are given in appendixes C 
and D. It can be shown from equations (5) and (6) that the boundary conditions for <p^ 
and $Q are 

Bcf) 




a 


dr 


(x,rb,co) = 0 


d(j) 


doj 


^(x,r,0) = -r 


^yg(x) £ r g yj(x) or r,^(x)^ 


J 


(46) 


and 


= cos a — 
8r I ' bM dx 


9$ 


8co 


(x,r,0) = 0 


(y 2 (x) g r s y^(x) or r^^(x)) 






(47) 


From equations (7), (8), (9), and (45), it follows that is governed by the equation 




A 


ax' 




a^0 


dX' 


+ 2M^“ sin^O! 


2 . 2 1 


dr dx dr p2 dco dx da) j 


(48) 
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where 


X = A 


(49) 


if 0^ is governed by equation (8) and 

X = A - ^1 - (50) 

if 0^ is governed by equation (9). It should be noted in passing that for near-sonic 
flow, equation (48) can be approximated as 


= (y + 1) sin^a ^ + 2 sin^oil 

^ u 3x 9x^ 


' 2 
30 3^0 



> 3r 3x 3r v,2 3 oj dx 3o;> 


30 3^0 


(51) 


It is assumed that in the region ^r < y 2 (^)] the potential 0^ in the governing 
equation (48) for can be approximated by the three-term expansion about the axis 

(r = 0) of the slender -wing potential for attached flow past a wing with no trailing edge. 
This approximation for 0^ is written as 




1 



COS 2 u> 


- r sin (jo 


(52) 


where the plus and minus signs apply above and beneath the wing, respectively. For wings 
with swept trailing edges, it is assumed that equation (52) applies ahead of the point where 
the trailing edge emerges from the body, and that the right side of equation (48) vanishes 
aft of this point. 


The particular solution to equation (51), and hence, for practical purposes, equa- 
tion (48), with (p^ approximated by equation (52) is 


^0 




^2 

r — - sin CO 
^2^ 



^ 2^2 



(53) 
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Since the complementary solution for is a source term (a term proportional to log^ r) 
and hence dominates the particular solution near the axis, it is concluded that there is a 
substantial region around the axis where the dominant part of is independent of w. 

In this report it is assumed that #Q(x,r,o>) can be approximated by a function 
which is independent of w. This function is defined as 

2 nC0=Tl/2 

^o(x,r) = - \ $Q(x,r,ci>) dw (54a) 


for points above the wing (a; > 0) and as 


(54b) 


for points below the wing (co < 0). If the computational plane is the wing plane (w = 0), 
values for (Pq{x,t) are calculated both above and below the wing. The governing equa- 
tion for <pQ(x,r) is 


A 


d^<P0 


dx^ 


1 

r ar 



+ sin a 



2M 


Sm Q! 




(55) 


The right side of this equation was obtained by integrating the right side of equation (48) 
with 0^ approximated by equation (52) over the regions above and beneath the wing. 
The equation for the function A used in equation (55) is 


A = 1 - - (y + 1)M^^ 


'dcpQ 


+ sin a 


u„(x) 


(56) 


The quantity u^(x) is the average velocity due to lift in the regions above and beneath 
the wing near the axis and is given by the equations 


u^(x) = ±y 2 ’ 


for attached flow, 


u (x) = ± 


1 + 


^2' + 


(57a) 


(57b) 
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for flow which separates at the leading edge with a vortex, and 


Uq,(x) = ±j 


, , ^ ^v)V 


yo + 


2y A T ' 

V V 

(l + 


(57c) 


for flow which is attached at the leading edge but which has a strake -generated vortex 
inboard of the wing tip. The plus and minus signs apply above and beneath the wing, 
respectively. Equations (57) are derived in appendix C. The boundary condition which 
is used at the axis is obtained from the first of equations (47) as 


Limit r 
r-0 \ 



= r^,(x) 


dr^,(x) 

dx 


(58) 


For free-stream Mach numbers below 1, the trailing-edge Kutta condition must 
be satisfied at least approximately in order for the solutions to be physically realistic. 
For a configuration with an unswept trailing edge located at x = x^-, this condition can be 
expressed as 


<;^?(x^+Ax,y,+0) - cp^x^+Ax,y,-0) = c/?^x^,y,+0) - (59) 

where Ax is the mesh spacing for the x-coordinate and z = ±0 designates the upper 
and lower surfaces of the wing and trailing vortex sheet. For shockless flow which is 
attached at the leading edge, this condition is enforced automatically if 0^ is obtained 
from the theory of Lawrence and Flax. If 0^ is obtained from slender-wing theory, 
the condition is met for shockless attached leading-edge flow by requiring that the deriv- 
ative Y2 approach zero at the trailing edge of the wing. It is shown in appendix C 
that if there is a vortex over the wing, the contribution of the vortex system to the deriv- 
ative ^0Q^’ on the inboard portion of the wing is of the form 

= (60a) 

if the flow is separated at the wing leading edge and 


A0 


a?,u 




(1 - \)V 

(1 - + X ^2 


(60b) 
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(60c) 



1 + Ty 




if the vortex is strake generated. Of course, equations (60) are not valid on the outboard 
portion of the wing near the vortex. If the computational plane is the wing plane where 
o) = 0, partial account of the effect of shock waves and leading- edge separation on the 
Kutta condition is taken by imposing the condition 


^0,u(V^^’^) - = ^0,u(^t’^) - ^0,f(^t’^) - " ^^«,u(^t) - 


(61) 


The subscripts u and H designate values on the upper and lower surfaces of the wing 
and vortex sheet. This additional adjustment for the Kutta condition is not made for 
computational planes other than the wing plane. 

Method of solution. - The method of successive line overrelaxation is used to solve 
equations (42) and (55). The procedure used here is very similar to that used by Bailey 
(ref. 28). However, the shock difference operator developed by Murman (ref. 29) is 
incorporated in the method since it tends to predict the jump condition at strong shock 
waves more accurately than other difference operators. The details of the numerical 
procedure are given in appendix F. The solution procedure used to find the vortex core 
strength and location is given in appendix G. 


Determination of Surface Pressures, Lift, 
and Pitching Moment 

Although the velocity potential is calculated in wind-oriented coordinates in this 
report, the surface pressure coefficients are calculated in body-oriented coordinates. 
The velocity potential is used in the form 


(p{x,r,e) = sin a cj)^{x,r,e) + (Pq{x,t) - | sin^a f'(x) (62) 

where the coordinates x, r, and 9 in this equation are body oriented. The angle 6 
is zero in the wing plane and positive above the wing. One of the forms discussed in 
appendixes C and D is used for the cross-flow potential cp^, and the solution determined 
numerically in terms of wind-oriented coordinates is used for the potential cpQ. The 
last term in equation (62) is a correction factor which accounts for the fact that the wing 
boundary condition for the lift potential 0^, is enforced, in the horizontal plane (o) = 0) 
when wind-oriented coordinates are used and in the wing plane {6 = 0) when body-oriented 
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coordinates are used. Let the body-oriented perturbation velocity components in the 
axial, radial, and cross-flow directions be u, v, and w, respectively. The pressure 
coefficient is written to second order in 6 sin a as 


= -2 u + sin a{v sin 0 + w cos 0^ - v^ - 


w 


(63) 


It was. shown by Munk (ref. 30) that the lift acting on that part of a slender configu- 
ration between the nose and the points x upstream of the trailing edge is 


L(x) = -cos a 



cp{x, y,z) dy + Trr^'^(x) sin 


“] 


(64) 


2 2 

where L is made nondimensional with the factor II where is the free- 

stream density. The path of integration C for the integral in equation (64) is the con- 
tour around the configuration in the cross-flow plane which intersects the axis at x, and 
the velocity potential cp is given by equation (62). Equation (64) can be expressed as 


L(x) = Lq(x) + Lj^(x) 


where 


Lq(x) = 2 cos a 


py=yo(x) ^ ^ 

\ A^o(^»y) ‘^y 

-y=0 


(65) 


and 


Lj^(x) = -cos a sin Q! ^ 4>^{x,y,z) dy + wr^^^(x) 


( 66 ) 


The quantity AcpQ(x,y) is the difference in the potential cp^ on the upper and lower 
surfaces of the configuration at the point x,y. The lift term Lq(x) in equation (66) 
is evaluated numerically in this report. 


It can be seen from the results of Wei, Levinsky, and Su (ref. 31) that for slender- 
wing flow past a configuration with a vortex over the wing and an unswept trailing edge, 
the lift term Lj^(x) can be written as 




Hx) 


y2(x) 


+ r, "‘(x) + 4v,(x) X.^(x) 


yg(x) 


y2(x) 


(67) 
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The functions and are associated with the vortex and are discussed in appen- 
dix C. As is pointed out in the appendix, these functions have different forms for leading- 
e(%e separation and for a vortex inboard of the leading edge. If there is no vortex over 
the wing, the last term on the right side of equation (67) vanishes. It should be noted that 
equation (67) applies from the nose of the configuration to the trailing edge of the wing. 

It has been shown by Mirels (ref. 32) that if the trailing edge as well as the leading 
edge of the wing is swept at x, the derivative dL^^dx for a configuration with a con- 
stant body radius can be written as 


dL^ 


(x) = 2ir cos a sin a 772 



( 68 ) 


where 772 given by equations (Cll) and K(k) and E(k) are complete elliptic inte- 
grals of the first and second kinds, respectively, with the modulus k given by equa- 
tion (C12). The fimction S(x) is discussed in appendix C. It can be seen that with the 
approximations given by equations (CIS) and (C16), equation (68) can be approximated as 


dLj 

dx 


(x) 


= 277 COS a sin a 



(69) 


where the derivative df/dx is given by equation (23). The last term in equation (69) is 
added so that this equation will be consistent with equation (67). The lift term Lj^(x) at 
points where there is a swept trailing edge must be determined by numerical integration 
of equation (69). 

Neither equation (67) nor (69) is valid at points aft of the trailing edge of the wing 
unless the body radius is constant behind the wing. Adams and Sears (ref. 33) derived an 
integral expression for the load on the aft end of a configuration with an unswept trailing 
edge and a body which closes. It is observed in reference 33 that because the flow in the 
aft cross-flow planes is essentially incompressible, the area of a stream tube with the 
outer and inner radii r^. and at the trailing edge is constant over the afterbody 

of the configuration. Consequently, the outer boundary of this stream tube is 

>■ = 

Aft of the trailing edge equation (66) for the lift term L^(x) can be written as 
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Lj^(x) = -cos a sin 


0Q,(x,y,z) dy + 


nr. 


“(X)] 


= cos a sin a 


,u}=2n 




+ 2 


pyt=y2(^t) 

V^b(^t) 




( 71 ) 


where is the difference in the linear lift potential <f>^ above and beneath the 

wing. The derivative dy^dy^ can be evaluated with equation- (70) . 

The pitching moment about the nose is calculated with the equation 

M = ^ dx (72) 

P Jx=0 dx 

TRANSONIC AREA RULE FOR LIFTING CONFIGURATIONS 

It is well known that the effect of thickness is to deflect streamlines outward from 
the configuration. It can be shown that for near-sonic flow, lift can also have this effect 
in addition to the usual downwash effect. The manner in which lift causes the outward 
deflection of streamlines is depicted schematically in figure 3. The cross-sectional 
area of stream tubes is minimum where the flow is sonic. The effect of lift is to increase 
the fluid speed in stream tubes above the wing and to decrease the speed in the tubes 


stream tube size 
near minimum in 
transonic free 
stream 




Velocity increases above wing; 
stream tube size increases 



Velocity decreases below wing; 
stream tube size increases 



Figure 3.- Outward displacement of streamlines by both lift and thickness effects. 
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beneath the wing. For configurations traveling at Mach numbers near 1, both the 
increase and decrease in fluid speed are deviations from near-sonic flow. Thus the 
cross sections of practically all the stream tubes about the body increase so that the 
streamlines are deflected outward more than they would be by thickness effects alone. 

It should be noted that this phenomenon does not occur for completely subsonic or com- 
pletely supersonic flow, where an increase in stream tube size on one side of the wing is 
compensated by a decrease on the other. 

An expression for the effective cross-sectional area of a lifting configuration 
traveling at transonic speed is given by equation (B32). The lift coefficient is defined in 
terms of the present notation as 


C 


L ^ 


2L 

^ref 


(73) 


where is the reference area. With equation (D23) and the fact that to lowest order 

the linear lift and the total lift L are the same, equation (B32) can be written to 
lowest order as 


/S dC \ 

Seff(x) = | C^{x) sin a (74) 

An expression for the coefficient B is given by equation (B33). The right side of this 
equation is at best only a weak function of x and can, for practical purposes, be con- 
sidered constant. However, a unique value for B cannot be determined with this equa- 
tion because the length scale £ cannot be determined with the theory of appendix B. A 
value for B can be obtained numerically with the present method by comparing the flow 
fields of lifting configurations and axisymmetric bodies with cross-sectional area distri- 
butions obtained from equation (74). 

It should be noted that the last term in equation (74) can be ignored for practical 
pvu'poses because both Cj^ and sin a tend to be much smaller than the gradient of Cj^ 
in the region where the effects of lift are large. It should also be noted that the results 
of references 13 and 14 indicate that there are higher order terms in the expression for 
the effective cross-sectional area which are proportional to the derivatives of the wing 
twist and camber. These terms are negligible unless the wing twist and camber are of 
the same order of magnitude as sin a. Since the wing is assumed to be flat in this 
report, the additional terms do not appear. 
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RESULTS 


In this section results are presented for axisym metric flow fields and flow fields 
about lifting configurations. In addition, the results of a numerical study of the transonic 
area rule for lifting configurations are presented. 

Axisymmetric Flow Fields 

Parabolic arc of revolution .- The present method calculates axisymmetric flow 
fields in essentially the same manner as the method of Bailey (ref. 28). The only major 
difference is that the present method employs the shock difference operator of Murman 
(ref. 29) and the method of reference 28 does not. In figure 4 results of the present 
method and that of Bailey are compared for flow past a parabolic arc of revolution with a 
fineness ratio of 10. The free-stream Mach number is 0.99, Results for the sonic-line 
and shock-wave locations and the surface pressure distributions are given. It is seen 
that the use of the shock difference operator influences the results for the location of the 
shock wave and the pressure in the region behind the shock. 


Shock-wave and sonic-line locations 



0 .2 .4 .6 .8 1.0 

x/C 

Figure 4.- Comparison of methods for parabolic arc of revolution with fineness ratio 

of 10. M^=0.99; a = 0°. 
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Parabolic arc of revolution with bump.- Equation (74) indicates that the flow fields 
about wing-body combinations traveling at near-sonic speeds at angle of attack are simi- 
lar to flow fields about bodies with bumps on them. Therefore, the present method should 
be capable of calculating flow about bodies with irregular thickness distributions. In fig- 
ure 5 the results of the present method are compared with the experimental results of 
Taylor (ref. 34) for a body with a bump on it. The basic body is a parabolic arc of revo- 
lution with a fineness ratio of 14 and the amount the bump, which extends from x = 0.393 
to X = 0.607, adds to the configuration radius is given by the equation 

’"/x - 0.393)1 
3 0.214 )] 

where is 20 percent of the basic body diameter. It should be noted that rg, 

the radius of the surface on which the boundary condition given by equation (58) is applied, 
was given the value of 0.03. This value is approximately the body radius where the bump 
begins and ends. The results for the pressure distribution differ slightly for other values 


Ar = Ar 


max 


sin‘ 



Figure 5.- Pressure distribution for bumpy body. = 0.975; a - 0°. 
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Ellipse of revolution.- The results of the present method for an elliptic body with a 
fineness ratio of 10 are compared with the results of the method of South and Jameson 
(ref. 35) in figures 6 and 7 for Mach numbers 0.99 and 0.997, respectively. In the present 
method the small -disturbance approximation is used, whereas in the method of refer- 
ence 35 the full potential equation is solved. It can be seen from figures 6 and 7 that, in 
general, the results of the two methods for the pressure distributions are in reasonable 
agreement. The results of the two methods differ in the region of the nose for = 0.99. 
This difference is probably due to the fact that the small-disturbance approximation, which 
tends to be inaccurate in stagnation regions, was used in the present method. It should be 
noted that the method of reference 35 does not employ the Murman shock difference oper- 
ator. Consequently, the results of the present method in the region of the shock are 
mathematically more realistic than those of the method of reference 35. It can be seen 
from the figures that the present method predicts a smaller supersonic region at large 
radial distances from the configuration than the method of reference 35. 

Parabolic arc of revolution in slotted and porous tunnels.- The results of the pres- 
ent method for transonic flow past a parabolic-arc body with a fineness ratio of 10 in a 
slotted- wall tunnel are shown in figure 8. The free- stream Mach number is 0.99. The 
body has a sting which extends downstream from x = 0.84. Comparisons are made with 
the results of the present method for flight in free air and the experimental results of 
Taylor and McDevitt (ref. 36). The experimental results were obtained in the Ames 


Shock -wave and sonic -line locations 
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-2h- 



Present method 

-- South and Jameson 

Body pressure distribution 
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Figure 6.- Comparison of methods for ellipse of revolution with fineness ratio 

of 10. M^ = 0.99; a = 0°. 
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Shock-wave and sonic-line locations 


Body pressure distribution 



x/a 

Figure 7.- Comparison of methods for ellipse of revolution with fineness ratio 

of 10. M^ = 0.997; a = 0°. 


Shock wave and sonic line Surface pressure distribution 
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Figure 8.- Parabolic arc of revolution in slotted-wall tunnel. M = 0.99' a = 0° 



14-foot wind tunnel. This tunnel has the same type of slots on all four walls, is essen- 
tially square, and has an openness ratio of 0.054 and a slot spacing of 24.77 cm (9.75 in.). 
The walls are actually neither slotted nor porous in the usual sense because the slots are 
subdivided with a corrv^ated strip. The theoretical length of the symmetric parabolic- 
arc body is 203.2 cm (80 in.). Results for two slotted-wall computations are shown. The 
dashed line depicts the results when the constant k in the slotted-wall boundary condition 
(eq. (29)) is evaluated with equation (30), which was developed by Baldwin, Turner, and 
Knechtel (ref. 26), and the characteristics of the Ames 14-foot tunnel. It is seen that the 
computation is in reasonable agreement with the experimental results except at the rear 
where the flow is probably separated. A comparison with the free -air results shows that 
the effect of the wall is to truncate the supersonic region and move the shock wave forward. 
It is known that the shock can be moved rearward by reducing the tunnel openness, and 
hence by increasing the value of k. It is seen that for the choice k = 1 the tunnel 
results and the free-air results are essentially coincident. Note that the shock location 
and the pressure level ahead of the shock are both in agreement. 

Results for flow in free air and a porous-wall wind tunnel are compared in figure 9. 

It is seen that the shock wave is located upstream of the free-air shock if the constant P 
in the porous-wall boundary condition (33) has a value of 0.5. This is the value of P 
which was recommended for the Ames 14-foot tunnel by Bailey (ref. 28). When P is 


Shock wave and sonic line 


Surface pressure distribution 




Figure 9.- Parabolic arc of revolution in porous-wall tunnel. = 0.99; a = 0°. 
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given the value 0.1, it is seen that the shock wave is located downstream of the free -air 
shock. However, the same minimum pressure level upstream of the shock, which differs 
from the free-air level, is obtained for both values of P. It appears from these calcu- 
lations that slotted-wall tunnels are more capable of reproducing free-air results than 
porous-wall tunnels. 


Flow Fields About Lifting Configurations 

Comparisons with linear subsonic and supersonic methods . - The results of the 
present method for incompressible attached flow past two wing -body combinations are 
compared with those of a standard linear panel method (ref. 37) in figure 10. The linear 
lift potential used in the present method is obtained from the theory of Lawrence and 
Flax (ref. 21). The body is an ogive-cylinder-ogive with each ogive being one quarter of 
the body length. A configuration which has a flat -plate delta wing with a 45° leading- 
edge sweep and a configuration which has a truncated flat -plate delta wing with the same 
sweep are treated. The results for the body pressure distribution in the wing plane, 
that is, the pressure distribution along the intersection of the body surface and the wing 
plane, are compared. The results of the two methods are in fair agreement. It can be 
seen that the panel-method results for the afterbody pressure depicted in figure 10(a) 
contain large oscillations and are probably in error. 

Present method 

Woodward-Carmichael 

Planform view of configuration Body pressure distribution in wing plane 




(a) Delta wing. 

Figure 10.- Comparison of results for body pressure distribution in wing plane. 

= 0; a = 4°. 
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Planform view of configuration Body pressure distribution in wing plane 
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(b) Cropped delta wing. 

Figure 10.- Concluded. 

In figure 11 the results of the two methods are compared for supersonic attached 
flow past two configurations. The linear lift potential which is used in the present 
method is obtained from a quasiconical theory. The shock-wave and sonic -line locations 
as predicted by the present method are shown on the left, and the surface pressure dis- 
tributions in the wing plane are compared on the right. Sonic-line and strong-shock- 
wave locations are determined from the longitudinal Mach number distributions. It is 
assumed that weak shock waves are located at points in regions of supersonic flow where 
the derivative cp" has a negative minimum. The free-stream Mach numbers for the 
delta-wing and the truncated delta-wing configurations are 1.1 and 1.4, respectively. For 
the lower free-stream Mach number there is a slightly detached bow shock, a tail shock 
which intersects the body just upstream of the rear apex, a weak shock in the vicinity of 
the wing, and small pockets of subsonic flow at the ends of the body. The results of the 
two methods for the surface pressure distribution are in good agreement in the vicinity 
of the wing. There is some disagreement on the forebody where the flow is transonic and 
hence where linear methods tend to be inaccurate. There is also some disagreement in 
the results for the afterbody pressure distribution. For the flow field with M^ = 1.4, 
the shock -wave pattern is qualitatively the same as that with M^ = 1.1, but there are no 
pockets of subsonic flow. The results of the two methods for the pressure distributions 
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(a) Delta wing. =1.1; 2°. 
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Shock-wave locations Body pressure distribution in wing plane 
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(b) Cropped delta wing. = 1.4; ot = 4*^. 

Figure 11.- Comparison of results for shock- wave and sonic-line locations and 
body pressure distribution in wing plane. 
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in the vicinity of the wing disagree somewhat, but the pressure difference across the 
wing is the same. The results for the afterbody pressure distributions are in agreement. 
The pressure increase on the rear part of the wing indicated in figure 11(b) is due to the 

first term on the right side of equation (55). This is the nonlinear ^ term as 

ax 

approximated in the present method at points inside the wing tip. As has been pointed out 
previously, this term is evaluated on the basis of slender-wing theory and, as a result, 
the Mach lines are assumed to be perpendicular to the free-stream direction. Of course, 
this approximation is poor for a Mach number of 1.4. 

Results obtained with the present method and the method of reference 37 for the 
chordwise pressure distributions at two span stations on the delta-wing configuration 
traveling at both subsonic and supersonic speeds are compared in figure 12. It is seen 
that the results of the two methods are in reasonable agreement. 

Transonic lift effects for subsonic free-stream speeds.- The effect of lift on the 
shock -wave and sonic -line locations in the wing plane of a wing-body combination traveling 
at a Mach number of 0.99 is shown in figure 13. The linear lift solution is obtained from 
slender-wing theory. The configuration is composed of a body with a fineness ratio of 12 
and the maximum-thickness point at x = 0.4 and a ''hyperbolic” wing. (See appendix A.) 
It is seen that for angles of attack of 6^ and less, the effects of thickness are dominant, 
and for angles of attack of 9^ and larger, the effects of lift are dominant. For the lift- 
dominated cases, the flow decelerates with a shock immediately ahead of the apparent 
bump due to lift, then overexpands past the bump, and finally decelerates with a shock 
near the trailing edge. 

In figure 14 the effect of lift on the sonic -line and shock-wave locations for flow 
past another wing-body combination traveling at the lower free-stream Mach number 
of 0.98 is shown. The configuration is composed of a circular-arc body with a fineness 
ratio of 10 and a flat -plate truncated delta wing with a rounded tip. As before, the linear 
solution is obtained from slender-wing theory. It is seen that the lift effect for this case 
is qualitatively the same as before, but that the angle -of -attack range for thickness- 
dominated flow is smaller. Lift-dominated flow exists for angles of attack of 6^ and 
larger. 

The pressure distributions along the intersection of the body surface and the wing 
plane of the configuration shown in figure 14 for a = 0^ and a - 9^ are compared in 
figure 15. For the larger angle of attack, the distributions both above and beneath the 
wing are shown. It is seen that the strength of the shock wave due to lift on the upper 
surface of the wing for a ^ 9^ is much larger than the thickness -induced shock for 
a = 0^. It is also seen that there is a small supersonic region and shock wave beneath 
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(a) M =0; 01 = 4°. 
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Present method 

Woodward-Carmichael 

Pressure distribution at y/b = 0.35 Pressure distribution at y/b = 0.65 



(b) = 1.1: a= 2°. 


Figure 12.- Comparison of results for chordwise pressure distributions at 
two span stations on a delta-wing configuration. 
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Figure 15.- Comparison of body pressure distributions in wing plane for two 

angles of attack. = 0.98. 

the wing. It should be noted that the results shown in this figure were calculated in the 
horizontal plane (co = 0^). 

As discussed in a previous section, the linear lift solution can be obtained from the 
theory of Lawrence and Flax if is less than 1, slender-wing theory if is 

near 1, and the modified theory of Carafoli et al. if is greater than 1. For high sub- 
sonic values of either slender-wing theory or the theory of Lawrence and Flax can 

be used since the latter reduces to the former as approaches 1. In order to demon- 

strate that both theories are applicable in this range, results obtained with them are com- 
pared in figure 16 for a free-stream Mach number of 0.98 and an angle of attack of 6^. It 
is seen that the results of the two theories are in good agreement except near the point 
where the wing emerges from the body. The differences in the pressure distributions in 
this region are due to the fact that some of the derivatives associated with the geometry 
are calculated analytically for slender-wing theory and numerically for the theory of 
Lawrence and Flax. Closer agreement can be obtained with a refined grid. 

In figure 17 the shock -wave and sonic-line locations for a? = 9^ calculated in the 
plane of symmetry {co = 90^, -90^) are compared with those calculated in the horizontal 
plane (co = 0^). The extent of the supersonic region in the upper part of the plane of sym- 
metry (o) = 90^) is about the same as that in the horizontal plane. However, it is seen 
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Present method with Lawrence -Flax theory 


Body pressure distribution 



Figure 16.- Comparison of wing -plane results obtained with slender-wing theory 
and Lawrence -Flax theory. = 0.98; a - 6°. 


Wing plane 


Vertical plane 



Figure 17.- Comparison of shock -wave 
and vertical plane. 


and sonic -line locations in wing plane 
= 0.98; a = 9°. 
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that the computation for oo = -90^ indicates that there is no supersonic bubble in the 
lower part of the symmetry plane so that the flow is entirely subsonic. This lack of 
supersonic flow is caused by the decelerating effect of the lift -dipole term on the flow 
beneath the wing. There is a difference in the results calculated for a; = 0^ and 
CO = -90^ since the former indicates the existence of the small region of supersonic flow 
beneath the wing. This difference is probably due to the simple flow model in the region 
r < y2(x) and indicates a limitation of the present method. 

Transonic lift effects for supersonic free-stream speeds .- In figure 18 the 
horizontal -plane results for the shock-wave and sonic -line locations for slightly super- 
sonic flow past a wing -body combination at angles of attack a - 0^ and a = 10^ are 
compared. The free-stream Mach number is 1.02. It is seen that two of the effects 

of increasing the angle of attack are to move the rear shock forward onto the trailing edge 
of the wing and to increase the size of the subsonic region behind the rear shock. 

The pressure distributions along the intersection of the body surface and. wing plane 
are compared in figure 19. In addition to moving the rear shock forward and increasing 
the size of the rear subsonic region, lift is seen to increase the strength of the shock 
above the wing. 

The results shown in figures 18 and 19 were obtained with a lift solution obtained 
from slender-wing theory. It is shown in figure 20 that virtually the same results can 
be obtained for =1.02 when the lift solution is obtained from the quasiconical theory 
of Carafoli et al. (ref. 22). This is to be expected since conical theory for supersonic 
flow reduces to slender-wing theory as the free-stream Mach number approaches 1. 

Lifting transonic tunnel flows. - The results for transonic flow past a lifting config- 
uration in free air and in a slotted tunnel are compared in figure 21. The Mach number 
and angle of attack are M^ = 0.98 and a = 12^, respectively. The slotted-tunnel com- 
putation is for a configuration with a span of 152.4 cm (60 in.) in an axisymmetric model 
of the Ames 14-foot transonic wind tunnel. The results show that the outer part of the 
supersonic region is truncated by the tunnel wall and that the presence of the wall causes 
the shock to be shifted forward. 

The effects of a solid wall on the shock -wave and sonic-line locations in the wing 
plane for slightly supersonic flow past a lifting wing-body combination are shown in fig- 
ure 22. The Mach number and angle of attack are M^ = 1,02 and a = 10^, respectively. 
The results for free air are shown on the left, and those for flow in the tunnel are shown 
at the right. It is seen that the bow shock is reflected at the tunnel wall and that the 
reflected shock intersects the aft shock in such a manner that the subsonic pocket at the 
rear of the configuration is enlarged considerably. The location of the reflected shock 
in the region of the intersection is not shown because of the difficulty in locating weak 
oblique shock waves. 
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Figure 18.- Comparison of shock-wave and sonic-line locations in wing plane 
for two angles of attack. = 1.02. 



x/H 

Figure 19.- Comparison of body pressure distributions in wing plane for 
two angles of attack. = 1.02, 
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Figure 20.- Comparison of wing-plane results obtained with slender-wing theory 
and quasiconical theory. = 1.02; a = 6^. 
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Figure 21.- Wing -plane results for free air and slotted-wall tunnel. 

= 0.98; a = 12°. 
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Free-air flow 
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Figure 22.- Effect of tunnel wall on shock-wave and sonic-line locations. 

M^=1.02; a = 10°. 

Leading-edge separation.- The results of the present method for the dependence of 
the lift coefficient on the angle of attack a for separated flow past a delta wing 

are compared with experiment and the results of other methods in figure 23. The com- 
parisons for low-speed flow are presented on the left side of the figure. The experi- 
mental data are those of Peckham (ref. 38). It is seen that the results of both Brown 
and Michael (ref. 20) and Smith (ref. 39) overpredicted the data, the latter being more 
accurate than the former, and the results of the suction analogy of Polhamus (ref. 40) are 
essentially in agreement with the data. It is also seen that the results of the present 
method overpredict the data to about the same degree as those of Smith. The suction 
analogy was not employed in the present method because the spanwise effects of vorticity 
cannot be determined with it in its present form. Results for supersonic flow with 

=1.96 are shown on the right side of the figure. The experimental data are those of 
Hill (ref. 41). It is seen that the trends for supersonic flow are the same as those for 
low-speed flow. The suction -analogy results for supersonic flow were obtained from 
reference 42. 

The effects of leading-edge separation on transonic flow past a lifting wing body 
are shown in figure 24. In this figure the results for attached and separated flow past a 
wing-body combination with = 0.98 and a = 6° are compared. On the left side of 
the figure it can be seen that the only effect of leading -edge separation on the shock -wave 
and sonic-line locations in the wing plane is the slight backward shift of the rear shock 
toward the trailing edge of the wing. The body surface pressure distributions in the wing 
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Figure 23.- Leading-edge separation for delta wing with aspect ratio of one. 
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Figure 24.- Effect of leading-edge separation on flow in wing plane. 

= 0.98; a = 6°. 
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plane are compared on the right side of the figure. It is seen that separation affects the 
pressure distribution both above and beneath the wing considerably. 

The coefficients for the total lift L = Lq + Lj and the linear lift Lj^ for attached 
and separated transonic flow past the configuration depicted in figure 24 are compared 
in figure 25. The Mach number is = 0.98 and the angle- of- attack range is 
0° < Q! < 9°. The nonlinear lift term Lq is due primarily to the presence of shock 
waves. It should be noted that the Brown and Michael model for leading-edge separation, 
which was used to obtain these results, overpredicts the effect of separation on the lift 
coefficient. It can be seen that the present results show a small increase in lift due to 
nonlinear effects. 

The results for flow past a configuration which has a strake -generated vortex over 
the wing are compared with the results for attached flow for a Mach number and angle of 
attack of M^ =0.98 and a = 6*^, respectively, in figure 26. The flow is separated on 
the strake and attached on the main part of the wing. The vortex which was generated by 
the strake is located over the inboard portion of the wing aft of the point where attached- 
leading -edge flow begins. It is seen that the shock -wave and sonic-line locations are 
essentially unaffected by the vortex and that the pressure distribution is affected appre- 
ciably only on the strake. 

The effects of leading -edge separation on the pressure distribution of the wing body 
depicted in figures 10(a) and 11(a) for Mach numbers appreciably different from 1 are 
shown in figure 27. The configuration is composed of a flat -plate delta wing and an 
ogive-cylinder-ogive body where the length of each ogive is one quarter of the body length. 
The angle of attack is 9°. Results for separated and attached flow are compared for 
Mach numbers 0 and 1.4. The pressure distributions which are shown are those along 
the intersection of the body surface and the wing plane. This distribution in the vicinity 
of the wing is the root-chord pressure distribution. It is seen that for incompressible 
flow the load at the wing root is larger for attached than for separated flow. However, the 
lift coefficients for attached and separated incompressible flow for 0 = 9'^ are 0.63 
and 1.14, respectively. The apparent discrepancy is due to the fact that the lift due to the 
leading-edge vortex is exerted mostly on the outboard portion of the wing. In figure 27(b) 
it is seen that the effect of the leading-edge vortex on the root-chord pressure distribution 
is also small for supersonic flow. However, the lift coefficients for attached and separated 
flow at these conditions are 0.80 and 1.02, respectively. 

Swept trailing edge .- The results for attached and separated flow past a configura- 
tion with a swept trailing edge are compared in figure 28. The free-stream Mach number 
and angle of attack are M^ = 0.98 and a = 4°, respectively. It is evident from the fig- 
ure that the rear shock wave intersects the body upstream of the swept trailing edge. The 
present method does not yield good results at points where a shock wave intersects a 
swept trailing edge. 
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Figure 25.- Total and linear lift for attached and separated flow 
about configuration depicted in figure 24. = 0.98. 
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Figure 26.- Effect of strake-generated vortex on flow in wing plane. 
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Figure 27.- Effect of leading-edge separation on body pressure distribution 
in wing plane of configuration shown in figure 10(a). a = 9^. 
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(b) M^=1.4. 
Figure 27.- Concluded. 
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Figure 28.- Comparison of wing-plane results for separated and attached flow past 
configuration with swept trailing edge. = 0.98; a = 4*^. 



Lifting Transonic Area Rule 

It has been shown that far from a lifting configuration traveling at near-sonic speed, 
the flow field is similar to that about an axisymmetric body with a cross-sectional area 
distribution which depends on both the lift and thickness distributions of the configuration. 
This modified area distribution is given by equation (74). As noted previously, the second 
term on the right side of this equation has a coefficient which must be determined exper- 
imentally or analytically. In figure 29 the shock-wave and sonic-line locations in the 
wing plane of a lifting configuration are compared with those of several axisymmetric 
bodies with cross-sectional area distributions obtained from equation (74) with several 
different values of the constant B. The Mach number and angle of attack are Moo = 0*98 
and a = 6°, respectively. It can be seen in figures 29 and 30 that the supersonic region 
for the axisymmetric configurations lies behind those for the lifting configuration. This 
phenomenon occurs because it is assumed in the present approximate method for lifting 
configurations that the lift solution is of the slender-wing type. Consequently, the char- 
acteristics of the lift solution are vertical. As a result, the present approximate method 
tends to locate those nonlinear lift effects which occur at large radial distances from the 
configuration slightly upstream of their actual location. It is seen that the radial extent 
of the supersonic region for the axisymmetric lift- equivalent body obtained with the value 
B = 3.8 is the same as that for the wing- plane results for the lifting configuration. 

The constant B should depend only on the configuration geometry for near-sonic 
flow; it should be independent of the angle of attack and Mach number. In figure 30, axi- 
symmetric results obtained with the value B = 3.8 and lifting-configuration results are 
compared for different Mach numbers and angles of attack. It is seen that reasonable 
agreement is obtained with this value of B at different angles of attack and excellent 
agreement is obtained at different Mach numbers. 

It can be shown that the coefficient B is dependent on the configuration geometry. 
The configuration depicted in figures 29 and 30 is composed of a flat-plate wing with the 
planform shown and a body of revolution with a parabolic-arc generator and a fineness 
ratio of 10. In figure 31 results are presented for a configuration composed of the same 
wing and a body with an elliptic generator and a fineness ratio of 10. The Mach number 
and angle of attack are M^ = 0.98 and a = 6°, respectively. It is seen that for this 
configuration the coefficient B has a value of 2.9. 
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Figure 29.- Determination of coefficient B in lifting transonic area rule. 

= 0.98; a = 6°. 



Figure 30.- Comparison of wing-plane results for lifting and equivalent axisymmetric 
flow past configuration with parabolic- arc body. 
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Figure 31.- Comparison of wing-plane results for lifting and equivalent axisymmetric 
flow past configuration with elliptic body. =0.98; a =6°. 

CONCLUDING REMARKS 

An approximate method for calculating transonic flow about lifting wing-body con- 
figurations has been developed. The nonlinear, three-dimensional physical problem is 
reduced to a two-variable computational problem which is solved with the method of 
relaxation. The Mach number range in which the method is applicable extends from zero 
to low supersonic and the angle -of -attack range extends to angles of the order of the con- 
figuration thickness-length ratio. The method accounts for shock waves, leading-edge 
separation, and wind-tunnel wall effects. A modified form of the transonic area rule 
which accounts for the effect of lift is developed. This effect is explained from simple 
physical considerations. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
November 26, 1975 
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APPENDIX A 


CONFIGURATION GEOMETRY 

The configurations which are treated in this report are composed of axisymmetric 
bodies and flat-plate wings. The body and wing geometry which were incorporated into 
the computer program (ref. 18) used to calculate the results presented in this report are 
discussed in this appendix. 


Body Geometry 

In this report the body shape rj^(x) is specified by equation (2), where t is the 
maximum body thickness and Fg(x) is the body thickness distribution. Four basic body 
shapes are treated. These are the ellipse with 

Fg(x) = 2\/x - x^ 


a general pointed body with 


Fg = C^l - x) - (1 - x)" 


where the maximum radius is located at 

1 

n-1 


-k 


and where 


C = 


or 


(1 - x^) - (l - x^.) 


F (x) = C(x - x") 


where the maximum radius occurs at 

1 

n-1 


X. 


■ir-k 
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and where 


c= 1 

Xc - Xc" 

a cone with 

Fg(x) = X 


and an ogive -cylinder -ogive with 

Fg(x) = 8x(l - 2x) 

(0 g X £ 0.25) 

F^(x) = 1 

(0.25 S X s 0.75) 

Fg(x) == 8(1 - x)(2x - 1) 

(0.75 S X S 1) 


The equations for the general pointed bodies were developed by McDevitt and Taylor 
(ref. 43). The basic body shapes are depicted in figure 32. As the figure indicates, a 
sting can be simulated at the rear. 


Wing Geometry 

The leading and trailing edges of the wing are specified by the function y 2 (x) and 
yj^(x), respectively. In this report the trailing edge is a straight line; it may be normal to 
the axis or swept. Two types of leading-edge shapes are considered. Both shapes sat- 
isfy the condition dygy^dx = 0 at the trailing edge so that the Kutta condition is satisfied. 

The first leading-edge shape, which is depicted in figure 33, is a portion of a hyper- 
bola. The leadir^ edge intersects the x-axis at Xj^, the wing ends at x^, and the wing 
semispan is b. The foci of the hyperbola lie on the line x = Xg, and the asymptotes 
of the hyperbola have the slopes Bj and -B^ and intersect at the point x = Xg, 
y = (1 + A)b. The equation for this leading edge is 


Y2(x) = b 


(1 + A) - 


Ia + (1 + 2A) 


/^3 - 

2 
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Figure 32.- Basic body shapes. 



Figure 33.- First leading- edge shape. 
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The second leading-edge shape, which is depicted in figure 34, is constructed of two 
or three straight-line segments which are connected with the transition function of the 


y 



Figure 34.- Second leading -edge shape. 


first kind of Grabau (ref. 44). Three straight-line segments are used if the wing has a 
strake; otherwise, two segments are used. As before, the wing ends at x = Xg. If the 
wing has a strake, the leading edge emerges from the body at x = Xj. For wings with no 
strake the point x = Xj^ is where the leading edge crosses the x-axis. The line seg- 
ments Yj and yjj intersect at the point 

Af - Ag ®2‘^1 " 

and the segments yjj and yjjj intersect at the point 
b - Ao 

""n = ^no = ^ 

The values of y 2 at Xj and Xjj are 

*'2(’‘n) “ '’C ■ ^ 2 ) 
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For a wing with a strake the equation for y„(x) is 


Y2(x) = Aj^ + Bj^x + 


(B2-Bi)(x-Xj) 

1 - expQkj^x - Xj)2 1 - exp0k2(x - Xjj)^ 


B2(x - xjj) 


where 


k - ^2 - ^1 
^ ^10^1 

k 

2 - bA2 


For a wing without a strake the equation for Y2(x) is 


Y2 (x) = A2 + B2X - 


- ^II) 


1 - expQk2(x - Xjj)] 


If the trailing edge is swept, the function yj(x) is 


yi(x) = b 


X - Xr 




where X2 is the point where the trailing edge intersects the x-axis. 
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REVIEW OF ANALYTICAL SOLUTION FOR TRANSONIC FLOW 
ABOUT LIFTING CONFIGURATIONS 

In this appendix the analyses of references 13 and 15, on which the approximations 
used in this report are based, are reviewed. These analyses, which use slender-wing 
and slender -body theory and which are performed with the method of matched asymptotic 
expansions, show the general form of the solution but do not completely determine it. 

It is well known that in the transonic flow regime disturbances propagate large dis- 
tances outward in directions perpendicular to the free-stream velocity vector. Conse- 
quently, the characteristic length scale in the radial direction must be large except in the 
immediate vicinity of the configuration, where radial disturbances are scaled by the char- 
acteristic radial dimension of the configuration. The regions close to the configuration 
and at large distances from the configuration are generally referred to as the inner and 
outer regions, respectively. For the present problem the radial length scales of the dis- 
turbances in the inner and outer region are written as I and i/v, respectively, where 
is a positive dimensionless parameter such that 

p «1 (Bl) 

and where £ is the longitudinal length scale of the disturbances. It should be noted that 
the same longitudinal length scale applies in both the inner and outer regions. The value 
of is of the order of the configuration length, but its precise value is not known. The 
radial variables in the outer and inner regions are related as 

r = pr (B2) 

In the outer region the perturbation velocity potential can be expressed as 

<;o(x,r,co) = e^$j^(x,r,a)) + ,w) + . . . (B3) 

where the gage functions satisfy the inequalities 

ej « 1 

^iH-1 « ^i 
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The relative magnitudes of v, and the quantity 1 - can be determined from 

a study of the governing equation (3). In the outer region this equation can be expressed 
as 


(l - M J) 






+ = e^{y + 1 ) 




ax 


8x^ 


(B4) 


Since equation (B4) must reduce to the standard linear forms for subsonic and supersonic 
flows, the two terms on the left must be of the same order of magnitude. Since the coef- 
ficient of the derivative must change sign in the transonic region, the term 

on the right must be of the same order of magnitude as the first term on the left. It 
follows that u, and 1 - may be related as 


1 - M 


V =ej_ = 


K 


(B5) 


where K is an order-one constant. Cheng and Hafez (ref. 15) argue that the similarity 
laws for flow about lifting configurations must reduce to those for flow about nonlifting 
configurations as the angle of attack is diminished. Consequently, i> must be propor- 
tional to 6. In this report it is assumed that 

1/ = 6 (B6) 


The relationship between the angle of attack a and the other parameters can be estab- 
lished only after the nature of the solution in the inner region is determined. 

It is shown in reference 13 that the perturbation velocity potential in the inner 
region can be expressed to second order in a and 6 as 


2 /l\ 2 

<p{x,r,oj) = 6 loggl-j<pQ^j(x) + 6 ^g(x,r,w) + sin a (b^{x,r,u)) 


. o 

+ 




(B7) 


The governing equations for the component potentials in this expansion are obtained fol- 
lowing substitution of the expansion into equation (2) and collection of terms of the same 
magnitude. It is found that all the potentials in the expansion (eq. (B7)) except sat- 

isfy the two-dimensional Laplace equation in the cross-flow plane. For example, 

. 0 
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The potential satisfies the Poisson equation 


2 90, 

V2^02 = (y + 1) — ^ + 2 


ax 


ax'^ 


'90J j 90J 9^0^ 

9r ax dv 1*2 0CO ax 9co i 


(B8) 


The potentials cp^ j and cp^ are the slender-body thickness potentials, and 0j 
is the slender-wing lift potential. It is well known that ^ can be written as 


n,l = 


(B9) 


and that for large values of r, (p^ and 0j can be written as 


<^5 = g 5 (x) + Fg(x) Fg’(x) logg r 

(BIO) 

0^=f(x)®'"^“ 

(Bll) 


where gg(x) is an unknown function of x, Fg(x) is the nondimensional thickness distri- 
bution of the equivalent axisymmetric body for the configuration, and f(x) is the axial 
distribution of the dipole strength for the lifting configuration. The prime denotes differ- 
entiation with respect to x. It is also well known that equation (BIO) for cp^ applies for 
all the way to the body surface if the body is axisymmetric, and that very near the config- 
uration surface 0j can be written as 

0j=m^(x,y)-z (B12) 


where the plus and minus signs apply on the leeward and windward sides of the configu- 
ration, respectively. The lift potential is discussed in detail in appendix C. 


Now consider the second-order potentials 2’ 
reference 13 that the pertinent solutions for 0^ ^ and 


02 , 1 ’ ^ 2 - 

02 1 are 


^2,2 ~ 

> 

*^ 2,1 ^ ^ 2 , 1 ^^^ 


It is shown in 


(B13) 


where the functions g 2 2 ^nd g 2 j are, at this point, unknown. It can be seen from the 
results of reference 13 that the particular solution for which is governed by equa- 
tion (B8), can be written as 
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(B14) 


where X = y + iz is the complex variable in the cross-flow plane and the asterisk 
denotes the complex conjugate. From equation ( 6 ) it can be shown that 02 must satisfy 
homogeneous Neumann boundary conditions at the wing surface. Consequently, the com- 
plementary solution for 02 can be written as 


2 2 4j Js=y^(x) 


p p 

909 ^ 809 ^ 

(x,S,+0) (x,s,-0) 


dZ 


dZ 


loggQs - X)(s - X*g ds (B15) 


where the function g 2 (x) is unknown. If the trailing edge of the wing is not swept, the 
lower limit of integration in equation (B15) is 0. It was pointed out by Cheng and Hafez 
(ref. 15) that the solution 02 ^ for attached flow exists only if the derivative is 

finite at the leading edge of the wing. Note that for this case it is necessary to bend the 
leading edge of the wing in order to enforce the condition on , In reference 13 it is 
shown that the solution obtained in reference 15 is valid to a first approximation for the 
case of leading-edge separation. As a result, it is possible to treat the problem of flow 
past flat-plate wings. 

It is shown in reference 13 that the sum of equations (B14) and (B15) can be written 
for large values of r as 




y + 


— f'(x) f"(x) ^2 logg^ r + cos 2 co) 


f(x) f’(x) 1 


~(1 - COS 2(d) 

r2 

(B16) 


where 


H(x) 


S=y2(x) 


m'(x,s) P.V. \ m'(x,t)/-^ + log |t - sA ds dt + C ' m(x,s) ds 

16 jt Js=-y2(x) ^t=-y2(x) - s e' >f 2 w 


(B17) 


The function m(x,y) in equation (B17) is written as 


m(x,y) = -| 


m^(x,y) - m 


.(x,y^ 


(B18) 


and the notation P.V. designates the principal value of the integral. It should be noted 
that with the aid of equation (Bll), equation (B 8 ) for 02 can be written for large values 
of r as 
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y + 1 f'f"/. 

2 p2 \ 


COS 


2o)) 



(B19) 


and that the particular -solution terms on the right side of equation (B16) can be obtained 
easily from equation (B19). 

The dependence of the angle of attack a and the gage functions €2 and eg in 
the outer expansion for cp on the equivalent thickness ratio 6 can now be determined 
by inspection. By assumption, the effects of lift and thickness in the outer region are 
comparable. Formally, this means that the leading term in the outer expansion must 
depend on both lift and thickness effects. The procedure is to write the inner expansion 
for cp for large values of the inner radial variable r in terms of the outer radial vari- 
able r and to require that the gage function of the largest lift term in the resulting 
expression be equal to = 6^. The gage functions of the second and third largest terms 
are then equated to eg and eg, respectively. From equations (B2), (B6), (B9), (BIO), 
(Bll), (B13), and (B16), it is seen that the inner expansion (eq. (B7)) can be written for 
r » 1 in terms of f as 


cp = 5^[gg(x) + Fg(x) Fg'(x) logg?^ + 6 sin o: f(x) 


+ sin^Q! logg^(| 


^ 2,2 


(x) + f'f”! 


+ sin^Q! loggfi 




+ H'(x) + 2i^f'f" log f 
2 ^ 


. 9 

+ sin^Q' 


gg(x) + H'(x) loggf 


+ 


Z + ^ f'f’f 2 log ^r + cos 
8 \ ^e 



+ . . . 


(B20) 


In reference 13 it is shown that if equation (B6) is to be valid, the term with the gage func- 
tion sin^Qf logg^^ij must vanish identically. It can be shown that if this were not the 

case, the terms in the inner expansion (eq. (B7)) of third and higher order in sin a when 

2 

written in terms of f for r » 1 would have gage functions larger than = 6 . As 
a result, it follows from equation (B20) that 

g 2 , 2 (x) = -^-^f’(x) f”(x) 
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sin a = 



(B21) 


where a is an order-one constant. The second and third largest gage functions in equa- 
tion (B20) are 


^2 = 


6 sin a = 



e 



a = 


a252 

‘-e(l) 


(B22a) 


(B22b) 


It can be seen from equations (B4) and (B5) that the governing equation for the 
outer potential is 

_ 0 

K— ^ + V ^4. = (y + 1)^^ (B23) 

8x2 ^ ^ 9x 3x2 


When the outer expansion (eq. (B3)) is substituted into the governing equation (3) for <p, 
it is found with the aid of equations (B5), (B6), and (B22) for the gage functions e e 2 > 
and eg that the outer potentials $2 and 4>g are governed by the equations 


and 


8^<^2 

8x2 


+ ’2^*2 


(y + 





^^*^3 - 2 

8x2 ^ 


(y + 





^8^ 
9x 3x2 


(B24) 


(B25) 


respectively. 

Equations (B23), (B24), and (B25) can be solved iteratively in the inner part of the 
outer region where r « 1. For example, equation (B23) is written as 


’2^«1 = 


(y + 1) 


84 >, 




8x 


- K 


8x 


67 



APPENDIX B 


The complementary solution to the left side of this equation is obtained. This solution is 
used to evaluate the terms on the right-hand side of the equation to render it a Poisson 
equation. Then the particular solution is obtained. This process can be continued until 
the desired number of terms has been found. It should be noted that the coefficients of 
the complementary solutions are arbitrary functions of x. These coefficients can be 
evaluated by matching the outer expansion with equation (B20). The outer potentials 
^3 3^re thus found to be 


= G^(x) + |Fg(x) Fg’(x) + f*(x) f"(x^log^ 


(B26) 


# = f(x) + . 


(B27) 


#3 = G 3 (x) + H'(x) loggf + >l^f'(x) f"(x)^2 log^^: 


(B28) 


The functions G^(x) and G 3 (x) in equations (B26) and (B28), respectively, cannot 
be determined with this method. Note that to lowest order the effect of lift is that of a 
source. The doublet effect appears in the second approximation. 

An expression for the effective cross-sectional area of a lifting wing -body config- 
uration can be obtained from the results of this appendix. The source strength of the 
potential for large values of the radial coordinate is proportional to the lengthwise deriv- 
ative of the cross-sectional area. From equations (B3), (B5), (B22), (B26), (B28), and 
others, it can be seen that the source strength is of the form 

S^F^F^’ + sin^a W + gin^a H' 

e e ^\5j 2 

Consequently, the effective cross-sectional area ^eff’ which is 2 77 times the integral 
with respect to x of the source strength, is of the form 

Seff=^rbV)+2^^1oge(| 

The first term on the right-hand side of equation (B29) is the actual cross-sectional 
area of the configuration. The transonic area rule as formulated by Whitcomb only 
involves this term. The second term is the dominant area-due -to -lift term which was 
first described in reference 12. The third term in equation (B29) is the higher order 


{sin a f'(x0 + 277 sin^a H(x) 


(B29) 
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area-due -to -lift term fovind by Cheng and Hafez (ref. 14). It can be seen from equa- 
tion (B17) that this term is composed of two integrals, the first a double integral and the 
second a single integral. With equations (B12), (B18), and (D25), the second integral can 
be approximated as 

ps=yo(x) 

\ m(x,s) ds « 2 tt f(x) (B30) 

^^s=-y2(x) 

From the form of equation (B30) it follows that the x-derivative of the function m(x,y) 
can be written as 

where mp is a slowly varying fxmction of x. For conical flow mp is independent 
of X. It follows from equations (B17), (B29), (B30), and (B31) that the expression for the 
effective cross-sectional area can be written as 


®eff^^^ " ^ ^ o ^ B(sin a f'(x0 + 2 sin^o! f(x) 


® (f ) * S C-1 "’p*’^’'’ (rh * ® l) 


A precise value for B cannot be obtained from equation (B33) since the length scale f 
is not known. It should be noted that the double integral is, at best, a slowly varying 
fimction of x. 
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LIFT POTENTIAL FOR NEAR-SONIC FLOW 


The slender-wing lift potential is used in this report for near-sonic flow. This 
potential is governed by the equation 


and satisfies the boundary conditions 


d(b 

B0 d(b 

— ^(x,r,0) = ^(x,r,7r) = -r 

de ’ dd 


> 

(V 2 2 r 2 yj or 


(Cl) 


where x, r, and 9 are body-oriented cylindrical polar coordinates. In this appendix 
this potential is derived for attached flow past configurations with swept leading edges 
and swept leading and trailing edges, and separated- leading- edge flow past configurations 
with swept leading edges. 


Attached Flow Past Configurations With 
Swept Leading Edges 

A cross section of the wing-body configuration is shown in figure 35. The solution 
for attached flow past this configuration was obtained by Spreiter (ref. 6) and Ward (ref. 7). 
Let the complex variable in the cross-flow plane be 

X = y + iz = rei® (C2) 
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Figure 35.- Cross section of configuration composed of 
circular body and flat- plate wing. 
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The complex potential for this problem is written as 


W(x,X) = i(x - Y(x,X)) 




The cross-flow perturbation velocity potential is 


0^(x,r,0) = Re[w(x,x3 


_ Ulb l2 20.V 2 


+ r4(l+-^ + 4r^4 2e+y2^ 1 + 




For large values of r, equation (C3) can be written as 

r r 

w,.,x)=.x 1- 


_ iX 1 V 2 fl + ""b \ 1 5^2 


1 - . 
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Thus, the perturbation velocity potential for large values of r can be written as 

sin 9 


0a(x,r,0) =f(x) 


(C6) 


where 


nx) = iy2^(i + ^ 

^ ' y2 


(C7) 


The quantity f(x) is the local dipole strength of the lifting configuration. 


Attached Flow Past Configurations With 
Swept Leading and Trailing Edges 

The derivative d<p^^dx for this problem has been derived by Mirels (ref. 32) and 
Mangier (ref. 24). It should be noted that in addition to the boundary conditions given by 
equations (Cl), the boundary condition 


80 

-5j2(K,y,0) = 0 



(C8) 


is also used. The solution is obtained in terms of the transformed complex variable 


Z = X - 


V 

X 


(C9) 


The 

that 


physical and transform planes are shown in figure 36. 


the derivative 


8% 

8x 8Z 


is written as 


It is shown in reference 32 


where 


8^W 
8x 8Z 


-i7?2^2* 




r 2 
^b 

^2 = y2 - T" 
^2 






(CIO) 


(Cll) 
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U-2ir^ 


(a) Physical plane. (b) Transform plane. 

Figure 36.- Cross-flow planes for attached flow past configurations with 
swept leading and trailing edges. 

The function S(x) is the function Mirels (ref. 32) developed for enforcing the Kutta con- 
dition at the trailing edge, and K and E are complete elliptic integrals of the first and 
second kinds with the modulus 



It should be noted that Mirels' S-function is identical to the function H used by Mangier 
(ref. 24). 

The complex potential W, and consequently the perturbation potential 0^, can be 
obtained by integrating equation (CIO) with respect to Z and x. However, as Mirels 
and Mangier point out, it is very difficult to perform these integrations and the lift can be 
obtained from equation (CIO) comparatively easily without first obtaining 0^. Conse- 
quently, 0^ and 90^y^9x are not determined at points near the configuration in this 
report. 

The derivative 90^^9x can be determined approximately in the outer part of the 
flow field. At large distances from the configuration equation (CIO) can be written as 


9^W 
9x 9Z 


- -i^2^2' ^ 


E(kf 

K(k)^ 


Consequently, the derivative 9W/9x can be approximated as 


9W 

9x 




E(k)1 

K(k)J 
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Therefore, for large values of r, the derivative 8(j)^^dx can be approximated as 


90 


9x 


“ = S(x) 


1 - 


E(k) 


K(kU 


^2 


\sine 


dx 


y2 


(C13) 


In this report two approximations are made which simplify equation (C13). First, 
it is assumed that the derivatives y^’ and yg’ satisfy the inequality 


yg’Cx) «yj'(x) 


(C14) 


Mirels (ref. 32) shows that as a result of inequality (C14), the function S(x) for a wing 
can be approximated as 


S(x) = 



(Cl 5) 


The second approximation is that of Mangier (ref. 24) for the function E(k)/K(k) which 
is written as 


||g=l- 0 . 85(1 -0-k2) (C16) 


With the approximations given by equations (Cl 5) and (C16), equation (C13) can be written 
as 


^ sin 6 
9x dx r 


(C17) 


where 


£=0.85y2 



| y2 - yi 

^2 + yi 


(CI8) 


Separated Flow Past Configurations With Swept Leading Edges 

In this report two types of flow fields with vortices are considered. The first is 
leading -edge separation, and the second is flow with inboard strake -generated vortexes. 
These flow fields are depicted in figure 37. 
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Cross section 


Cross section 







Top view 



(a) Leading- edge separation. 

Figure 37.- Flow fi 


Top view 



s with vortexes. 


The problem of leading-edge separation was first treated successfully by Brown 
and Michael (ref. 20) for the case of flow past a delta wing. Although several other more 
complicated leading-edge-separation m.odels have been developed and could have been 
used, the model of reference 20 is employed here because of its simplicity. 

The complex potential for the Brown and Michael model for separated-leading-edge 
flow past a flat-plate-wing, axisymmetric-body configuration is 


W(x,X) = i / X - Y(x,X) - y(x) log^ 


Y(x,X) - Yv(x) 


Y(x,X) + Y *(x) 


(C19) 


where X and Y are given by equations (C2) and (C4), respectively, and where 


Yv(x) = 

is the vortex location in the Y-plane. The potential in equation (C19) satisfies the bound- 
ary conditions given by equations (Cl). Brown and Michael show that in addition, a Kutta 
condition must be satisfied at the leading edge and the force on the system composed of 
the vortex core and vorticity feeding sheet, which extends from the leading edge to the 
core, must be zero. These conditions are expressed as 


rb^(x) 



\{x) + iTv(x) 


(C20) 
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(C21) 


and 


dyy 


dx 


(Xv - y2) 


= 


Limit 

X-^Xy 



+ 


1 \ 
i X* _ x^* ) 


dx 


(C22) 


respectively. For the problem of flow with a strake -generated vortex, the Kutta condi- 
tion, equation (C21), is not enforced and the vortex strength becomes constant so that the 
left side of equation (C22) vanishes. The solution procedure is given in appendix G. 


The problem of separated flow past a delta wing can be used to gage the effect of the 
vortex in the outer part of the flow. It is shown in reference 20 that to lowest order in 
sin a the location of the vortex in the X-plane is 


Vy 5^2 


2lyo 


2/3'' 




z, = - X sin a 

V 4 


(C23) 


It follows that in the outer part of the flow where 


X|» Y2 


the perturbation velocity potential is 



2 sin a 
r 


1 + 2 



(C24) 


Consequently, for small angles of attack the second term in equation (C24), which is due 
to the vortex, could be ignored in the outer part of the flow. However, this term is not 
ignored in this report because its inclusion presents no particular difficulty. 

The influence of the vortex on the flow far from a general configuration can be 
determined from equation (C19). For large values of r and hence large magnitudes 
of Y, the logarithm in equation (Cl 9) can be approximated as 
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loge 


'y(x,X) - Yy(x)' 
Y(x,X) + Y/(x) 


Y^(x) + Y/(x) 
Y(x,X) 


2Ay(x) 


72 (x) 


V(x)i 

7o(x) 




Consequently, the potential <p^ at large distances from the configuration is given by 
equation (C6) where 


f(x) = I 72 





(C25) 


If the flow is separated at the wing leading edge, the boundary condition given by equa- 
tion (C21) applies. It can be shown that this boundary condition can be written in the 
form 



(C26) 


If the flow is not separated at the wing leading edge but a strake -gene rated vortex is 
present, does not vary with x aft of the strake and has the value given by equa- 
tion (C26) at the time that leading-edge separation ceased. 

In order to account for nonlinear behavior at points where r < y 2 (x), it is necessary 
to account for the effect of the derivative on the coefficient of the 9^^Qy4x2 

term in equation (55). In this report the value u^, of this derivative which applies in 
the region near the axis is used to approximate the values at all points inside the wing 
tip. From equation (C19) it can be shown that if r^ « y 2 , the potential 0^ near the 
axis is, approximately. 


0Q, = ±72 (x) + y.^,(x) 


tan"^ 




Consequently, the derivative u^ 


9x 


in the region near the axis is 


u 


Of 


±72* (x) 




, 2yy(^l - Ty ) X ^' ^Vv V 

(±1 - (±1 - 


(C27) 
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The function Y^(x) given in equation (C20) is related to the vortex core location 
Xy. = + iZy by the equation 


Y^(x) = 


(C28) 


It can be shown from equations (C20) and (C28) that for 
and are 



V. 


r^ « yg, the equations for 



(C29a) 



If 1 _ y^y^y 2 « 1 3,nd \jy2 case for separated-leading-edge flow past 

a wing at a small angle of attack, it follows from equations (C29) that 


(C30) 

With the first of equations (C23) it can be shown that for the particular case of separated 
flow past a delta wing, equations (C29) can be written as 



If 1 - ^v/y2 « 1> as is the case for a strake- generated vortex over 

the wing, the values of Xy and are of the form 


\ «1 


Ty « 1 
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\ 


V 





« 1 


(C31a) 







(C31b) 


With equation (C26) and inequalities (C30) it can be shown that if the flow at the 
leading edge is separated, equation (C27) can be approximated as 




^2’ 


1 + 


(C32) 


With equations (C31) it can be shown that if there is a strake -generated vortex over the 
wing, equation (C27) can be approximated as 


U^ = ±! 


^ A - \2 2 

(1 + 


2y A T ^ 

V V 


(l ? ^v) 


2 2 


(C33) 


In order to satisfy the Kutta condition to lowest order, it is necessary to correct 
for the influence of the vortex on the flow at the trailing edge. For all the slender-wing 
solutions used in this report, the attached-flow portion of the trailing-edge Kutta condition 
is satisfied by the condition that y 2 * vanish at the trailing edge. Consequently, it fol- 
lows from equation (C32) that if the flow at the leading edge is separated, the influence of 
the vortex on the flow on the inboard portion of the wing is 




= ±- 


1 + T, 


■(Vv)' 


(C34) 


It can be seen from equation (C33) that if there is a strake -generated vortex over the 
wing, the influence of the vortex on the flow on the inboard portion of the wing is 
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2y. 


2 2 


This equation can be simplified as 


|(±1 - 


A0' =2y ~ 


A0’ 


2y,; 


1 + T„ \ 




1 + T,, 


(C35a) 


(C35b) 
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LIFT POTENTIAL FOR SUBSONIC AND SUPERSONIC FLOW 

In this report approximate methods are used to determine the lift potential <p^ 
for subsonic and supersonic flow. These methods are described in this appendix. In 
addition, a generalization of the Brown and Michael model to subsonic and supersonic 
flow is given. 


Solution for Subsonic Flow 


The approximate theory of Lawrence and Flax (ref. 21) for subsonic flow past wing 
bodies is used in this report. This theory is a generalization of the theory of Lawrence 
(ref. 25) for wings alone. The theory of references 21 and 25 is restricted to wings with 
unswept trailing edges. 


The transformation given in equation (C9) is used to relate the flow past a wing 
body to that past a wing alone. These flow fields are depicted in figure 38. The quan- 
tity 7 } 2 is related to Y 2 and r^ by the second of equations (Cll). 


Z 



1 Z -plane 

2ir^ 


^2 




U-2ir^ 

(a) Wing-body plane. (b) Wing plane. 

Figure 38.- Wing -body and wing cross-flow planes. 


Lawrence employs the equation 


. „ ps=x. pt=T72(s) u (s,t) 

w(x,r?) = \ M ^ dsdt-y— - 

2n dj] ds=Xn •-^t=-no(s) H ~ t 


^2^s) 




X - s 


(Dl) 


to determine the downwash on a wing. The foremost points on the wing and trailing edge 
of the wing are at Xq and x^., respectively, and u„ is the x component of pertur- 
bation velocity on the upper surface. Equation (Dl) is an exact integral solution to 


I 
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equation (9) which was developed by Reissner (ref. 45). Lawrence evaluates this equation 
by using the crucial approximation 


^(x - s)2 + (l - mJ) (y - t)2 « I |x - s| + |(x - s)^ + (l - mJ) 




(D2) 


As noted in reference 25, this' approximation was chosen because equation (Dl) reduces to 
the limits of slender -wing theory and two-dimensional theory when rjg approaches zero 
and infinity, respectively, when the approximation is used. The boundary condition on the 
wing is 


w(x,y) = -a 


(D3) 


With equations (D2) and (D3), equation (Dl) can be written as 


S-Xq 


1 - s| +\/(x - s)^ + (l - 


2(x - s) 


9 f>t=7,2(s) Ua(s,t) 

»2(s) 


dV Jt=- 


dt ds (D4) 


Equation (D4) is used in the integral form 


P??=t? 2 (x) I — 27: r j ps=x 


|x - s| +|/(x - s)^ + (l - m„^)i)2^(x) 

1 + 1 > i—i! 


2(x - s) 


X 


P?7=?72(x) j 5 ;■ 3 /-t=7?2(s) u,,(s,t) 

, \/'?2 (D5) 

^ 97 Jt=-r7,(s) 


It is shown by Lamb (ref. 46) that any two potentials and ip 2 which satisfy 
the two-dimensional Laplace equation 

VjVi - 0 


V 2 V 2 = 0 


also satisfy the integral equation 



(D6) 
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where S is a closed path in the plane and n is the coordinate normal to the curve S. 
For this problem the path S is chosen along the surface of the wing so that the normal 
coordinate is simply Let the potential surface and its normal derivative 

be 


> 



(D7) 


It can be shown with slender-wing theory that and its normal derivative can be 

related as 


Jt=-7]2(x) 7 - t 


dn 


(D8) 


For this problem let 


= u(x,? 7 ) (D9) 

With the choices given in equations (D7) and (D9) for 1^2 and respectively, and 
with the relation given by equation (D8), it can be shown that integral equation (D6) for 
this problem can be written as 


p77=^72(x) 

'J77=-??2(x) 


/ 


2 o’ a rt=7o(®) u (s,t) rri=ri 2 i^) 

- 7 ^ d77 — \ — dt ds = 77 \ u„(s,t 7 ) d^ 


9 p77=7?2W 

= 77 sin q; — \ <pAs,v) dr? 

(DIO) 


where is the upper-surface value of a lift potential which, in general, is related to 

the potential y by the equation 


^a,u ~ r 


sin or 


(Dll) 
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Let the functions k(x) and g(x) be defined as 


k(x) = y? 72 ^(x) - 77^ d 77 = I 7 ? 2 ^(x) 

*^r7=-T72(x) '' 


‘^77=-772(x) ^ ^n=-r 


0u(x,7?) - 0 £(x,t]) d?7 




•^r7=-772(x) 

With equations (DIO), (D12), and (D13), equation (D5) can be written as 


k(x) = I 


2 g(x) 


ns=x^\l{x - s)2 + (1 - uj)r]^^{x) 

^ 


(s) ds 


Let the variables a and t and the function ^(a) be defined by the equations 


X = 


r ^ 


Ft ^0 (^t ■ ^0) ^ 


L ^ 

> 

X^ + Xq + ^x^ - Xq^ COS T 

J 


and 


3(<^) = 


Xt -Xq 

With equations (D15) and (D16), equation (D14) can be written as 


4k(cr) = 2g(a) + g(CT) + ^(a) ^ ■ + f H(ct,t) g'Cr) dr 

'^T=0 cos T - COS a 'Jj=0 


(D12) 


(D13) 


(D14) 


(D15) 


(D16) 


(D17) 


where 
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Equation (D17) is an integral equation for the function g(a). This equation is solved 
numerically by colocation. The procedure used is that outlined in reference 25. The 
form assumed for g(a) is 


N-1 

g(a) = (tt - a)Ai + ^ (^r-i “ Aj,+ i 

r=l 


sin (rg) 
r 


(D19) 


where 


■^0 “ -^N-1 " -^N ■ ° 

It will be shown that as a result of the form of equation (D19), the trailing-edge Kutta 
condition is satisfied to lowest order. 

The function g(x) is related to the load on the configimation between the apex and 
the station x. The linear lift is given by equation (66). With equation (Dll) it can be 
shown that the contour integral in equation (66) can be written as 


^ dy = 0(x,y,z) dy - 277rjj^(x) 


(D20) 


Let the contour C lie along the body and wing surface. It is shown in reference 21 that 
the integral along the body surface is related to that along the wing as 


pO)=277 ^y=y r 

\ dy = 2 \ 0jx,y) - <pAx,y) 

'^o)=0 ' ' -'v=r^(x) ^ 


y=Y2(x) 

y=r^,(x) 


r,.^(x) 


dy 


Consequently, the contour integral on the right side of equation (D20) can be written as 

dy = 2g(x) (D21) 


r / 1 

r, 2(x)' 

0u(x,y) - 0^(x,y) 

y^ 


Thus from equations (66), (D20), and (D21), it is seen that the linear lift Lj^(x) is related 
to g(x) as 


Lj^(x) = cos OL sin a. 


2g(x) + 7rr 


b^(x)] 


(D22) 


With equations (11), (12), and (66), a similar relationship between Lj(x) and the 
function f(x) can be established. This relationship is 
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L]^(x) = 77 COS Of sin a [2f(x) - 

From equations (D22) and (D23) it follows that f(x) and g(x) are related as 


(D23) 


g(x) = 77 


f(x) - r 


b^w] 


(D24) 


and from equations (Dll), (D13), and (D24), it follows that f(x) is related to the integral 
of wing surface potentials <f)^ and 0^ ^ as 




(D25) 


The method of Lawrence does not provide a means for determining the pressure at 
isolated points. In this report approximate values for the pressure are obtained by using 
the assumption that the spanwise pressure distribution is functionally the same as that 
obtained with slender-wing theory. The linear lift on the configuration can be written as 


H s=x p77=77p(s) 
\ 

3=Xq -’7J=_77 (s) 


- Cp 




ds d?7 + 77 sin a rjj^(x) 


(D26) 


A comparison of equations (D22) and (D26) shows that 
« pS=X p7?=77p(s) 

g(x) = ^ \ \ C Js,ri) ds dr? 

sin a '^s=Xq '^ri=-ri 2 {s) 


(D27) 


The spanwise integral of the pressure coefficient is obtained from equation (D27) as 


pr?=772(x) 

'^77=-tj2(x) 


Cp d, = I sin c g 


(D28) 


The pressure coefficient C . on the wing alone is expressed in a form proportional to 

Pj ^ 

that of slender-wing theory as 


Cp,,(ii,.|) = 


Rq sin a 


(D29) 
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where Rq is to be determined. When equation (D29) is substituted into equation (D28), 
the expression for Rq is foxmd to be 

The pressure coefficients on the upper and lower surfaces of the wing of the wing-body 
configuration are approximated as 



(x,y) 




dg/dx 



(D31) 


The expression for the derivative dg/dx can be obtained from equation (D19) and 
the first of equations (D15) as 


N-2 

dx X. - Xf. r 

t 0 r=l 


sin (ra) 


It can be seen that dg/dx vanishes at the trailing edge where a = 0. Thus, the approx- 
imate expression for the wing-surface pressure coefficient given in equation (D31) van- 
ishes at the trailing edge so that the Kutta condition is satisfied. 


Solution for Supersonic Flow 

The quasiconical theory of Carafoli et al. (ref. 22) is used to evaluate the linear 
lift potential for supersonic flow. As in the subsonic case, the wing body is mapped to a 
wing alone with the transformation given in equation (C9). This mapping is depicted in 
figure 38. 

The form which is assumed for the potential is given by equation (26). The 

form of the constant A used by Carafoli et al. is 


bf +b2 

^2{H) 

i 

+^2'(^2) 


E(k) 


(D32) 
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where the distances and bg and the positions and X 2 are depicted in fig- 
ure 39 and where E(k) is the complete elliptic integri of the second kind with the 
modulus 


k = 





y y 



(a) Geometry of Carafoli et al, (b) Present geometry. 

Figure 39.- Wing parameters for supersonic solution for (j>^. 


It can be seen that the value of A given by equation (D32) becomes small as ^2'(^l) 
becomes small so that the potential <p^ becomes large. In order to avoid this difficulty, 
a form for A given by the equation 


where 


A = 


2ci 


E(k) 
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is used in this report. The distances Cj and C 2 are depicted in figure 39. It can be 
seen from the equation for k that a real value for k exists if the wing leading edge is 
sonic or subsonic. Consequently, the present method is limited to flow about configura- 
tions with sonic or subsonic leading edges as determined by simple sweep theory. 

Brown and Michael Model for Subsonic and Supersonic Flow 

The leading-edge separation models of both Brown and Michael (ref. 20) and Smith 
(ref. 39) are for lift potentials governed by equation (8), the Laplace equation in the cross- 
flow plane. In this section the Brown and Michael model is generalized for lift potentials 
governed by equation (9) for subsonic and supersonic flow. It should be noted that Sacks, 
Limdberg, and Hanson (ref. 47) have generalized both the Brown and Michael model and 
the Smith model for subsonic flow. 

The potential <p^ is written as 

<^0! = <^a + ^v (1^33) 

where (p^ and are the attached flow and the vortex potentials. It is assumed that 

these potentials are governed by the equations 



and 


= 0 (D35) 

The form of equation (D35) is based on the fact that the cross-flow derivatives of 
are dominant. The complex potential which satisfies equation (D35) is 


Wy = -i7v(x) logg 


Y(x,X) - Y^(x) 
Y(x,X) + Y/(x) 


(D36) 


where X, Y, and are given by equations (C2), (04), and (020), respectively, and 
where y^(x) is the strength of the vortex core. 

The complex velocity associated with the attached-flow portion of the linear lift 
potential is 
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8y 8z 


(D37) 


The boundary conditions which are used to determine the strength and location of the vor- 
tex core for a separated-leading-edge flow are the leading-edge Kutta condition and the 
no-force condition on the system composed of the vortex core and the vorticity feeding 
sheet. These conditions are expressed as 



+ 


dX / 

''x=Y2 


= 0 


(D38) 


and 


dx 


(Xy - V2) ' 


Limit (v„* 


+ 


/V 1 \ 

dX* ^ i X* - X^*j 


<■ 

dx 


(D39) 


For the problem of a strake-generated vortex, equation (D38) is not enforced. 

The same general procedure is used to evaluate the derivatives and 

for subsonic and supersonic flow. It is assumed that the variation of 4>^ in the 
cross-flow plane is the same as that of slender-wing theory. Consequently, the poten- 
tial ip^ is written as 


Re 



(D40) 


where X and Y are given by equations (C2) and (C4), respectively, and where Rj is 
a function to be determined. With equations (D24), (D25), and (D40), it can be shown that 
Rj^ has the value 


7I772'^(x) 

Thus, for subsonic and supersonic flow the complex velocity is 


(D41) 


V -i 

^ 7rr?22(x) 


1 - ^ (x,X) 
ax ’ 


The solution procedure used to find the vortex core strength and location is given 
in appendix G. 
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APPROXIMATE FORM OF FAR- FIELD LIFT POTENTIAL 
FOR SUPERSONIC FLOW 


Heaslet and Lomax (ref. 48) show that the linear lift potential for supersonic flow 
past a wing can be written in integral form as 


^(x,y,z) = 


sin a 

2-n 


ps=x-/3r 

's=0 


(x - s) ds 


t=Y2(x) 

YgW (y - t) 


c 


2 


^u^(s,t) 
~R 


dt 


where Au^ is the jump in the x-derivative of the lift potential across the wing 

and where 


R 


= j(x - s)2 - 


(Y - t)2 + z2 




(El) 


For large values of r the potential c;^)j^(x,r) in the outer expansion (eq. (10)) can be 
written as 


^ (x,r) = - = i - 

1 smasmu, r ^(x . 


ds 


where 


f'(x) = ^ 

27T »Jt= 


Au (x,t) dt 

Y2W 


If f'(0) = 0, equation (E2) can be integrated by parts to yield 


^j(x,r) = p r ^ ^(x - s)^ - /3^r^ f”(s) 
^ r -^s=0 ^ 


ds 


If the flow is attached at the leading edge, f(x) is given by the equation 

fw = I 


(E2) 


(E3) 


(E4) 


91 


APPENDIX E 


From the results of appendix C it is seen that equation (E4) is also a reasonable approxi- 
mation for separated-leading-edge flow for the ai^le-of -attack range considered in this 
report. For wii^s for which the derivative 72 * (x) is either constant or a slowly varying 
function of x, equation (E3) can be integrated to yield 


<pj(x,r) = 







(E5) 


Since 72 ' can be approximated as 


yj'W 


y2W _ y2<* - M 

X X - )3r 


equation (E5) can be written in the form 


<;i’l(x,r) 


^(x - /3r) 



loge 
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NUMERICAL METHOD 


In this report, the method of successive line overrelaxation is used to solve equa- 
tions (42) and (55). Equation (42), which governs the flow at the points r > Y 2 (x), is used 
in the form shown at the points r > ^^^(x), where rg^(x) is slightly larger than Y 2 (x). 
Equation (55), which governs the flow at the points r < Y 2 (x), is used in the form given. 

In the region Y 2 (x) < r < i‘^(x), it is assumed that the outer expansion (eq. (10)) is appli- 
cable and that the terms proportional to sin co and cos 2co are governed by equations 
of the form of equation (8) if (p^ is obtained from slender-wing theory and by equations 
of the form of equation (9) if (p^ is obtained from the theory of Lawrence and Flax or 
from quasiconical theory. Consequently, the equation which governs (Pq for 
Y 2 (x) < r < r^(x) is 


^‘^0 1 d f 

2 r 3r \ dr J 



. 9 

sin^o? 


y + 1 fT^ 

2 r2 



sin a f- sin co h- 

r \ 2 r,2 ^4 


cos 2u) 


(FI) 


where 


A = 1 - mJ - (r + 1)M J 


cpf^' + sin a — sin o) + sin^af ^ ^ ^ ^ ^ 
^0 r I 2 r4 


cos 2ci) 


The function X is related to A by equation (49) or (50). 
Equations (42), (FI), and (55) can be written as 


A 


, 


dx^ 


ll^/r M) = F(x,r,o)) 
r ar \ 9r / ' ’ ’ ' 


(F2) 


where 


A = 1 - M J - (y + 1)M^ 


dip 

dx 


+ G(x 




(F3) 
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The function \}/ is equal to <p for r > and to cpQ for r S r^^, and the functions 
F and G are the appropriate known analytical expressions. The grid system used is 
virtually the same as that used by Bailey (ref. 28). In the radial direction the 
transformation 


r - ^0 + ^ 


arj 


is used so that r varies from rQ to as y varies from 0 to 1. The effect of 

the transformation is to concentrate points near the axis. The constants b and a 
must be appropriately chosen to give the desired point distribution. Note that if b = a, 
the outer boundary is at infinity. In the vicinity of the body a uniform grid is used for 
the x-coordinate. For = 1, the x-grid is stretched downstream and upstream of this 
region in the manner used in reference 28. This stretching continues the x-grid to the 
finite values and For > 1 the x-grid spacing is uniform throughout. 

For = 1, the boundary condition which is applied at x . and x „„ and at 
the outer radial boundary if that boundary is at infinity is cp = 0. For > 1, this con- 
dition is applied only at since no downstream boundary condition is required and 

since the outer boundary is located at a finite value of r if the free stream is super- 
sonic. The boundary condition which is applied at the inner boundary is given by equa- 
tion (58), and the one which is used at the outer boundary if that boundary is to represent 
a tunnel wall is given by equation (28) if the wall is solid, by equation (32) if it is slotted, 
and by equation (33) if it is porous. 

The partial derivatives with respect to x are approximated with central finite- 
difference expressions if the local Mach number M < 1 and backward expressions if 
M > 1. Let the indices for the x- and r- coordinates be j and k, respectively, and 
let 


- ’'j-2 

AXj = Xj - Xj.i I 

AXj = Xj^i - Xj 

J 

The central-difference expressions are 


(F4) 



Axo 


AXo (AXo + Ax, 


AXr 


AXri 


^ j '*■ Axg 


Ax, 


Ax2(ax2 +AX3) 




(F5) 
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V ■ TT 7Z~ t + T V (F6) 


^9x^/c AXg ^AXg + AXg^ j+ljk Ax 2 Axg AXg ^AXg + AXg^ 

and the backward expressions are 




8^iA ^ 2 _ 2 , 2 

AXo ^AXj + Ax 2^ ^^9! AXi ^AXi + Ax^^ 


■** ^^^2/ — ^ AXj + ^^2^ 


At strong shock waves the nonlinear term A^— ^ in equation (F2) is replaced by the 

dX^ 

Murman shock operator (ref. 29) which is written as 


/a®M - a r^A +A r n 


The partial derivatives with respect to the radial coordinate are written as 


r ar(^ '^k’^j,k-l " ®k'^j,k ■*■ ^k’^j,k+l 


(F7) 


where 


■teJ-;---? 


^ b^(A?7) 


7k + (b - a^kj^o 


- 2a 


Bk 


b2(A?7)2 


_ - ^k) ( An 

b2(Ar?)2 


7v + 


(b - ank)^ 


2a 


(F8a) 


(F8b) 


(F8c) 
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except at points on the boundaries. At the inner boundary the radial difference expres- 
sion is written as 



and at the outer boundary as 


1 

r dr 





km- I 


where the derivative 



is obtained from equation (28), (32), or (33). 


When the computational plane is the horizontal plane (oj = 0), the velocity potential 
is calculated both above and beneath the wing for r < y 2 (x). Modifications to the finite - 
difference formulas need to be made at several points neighboring the wing. If the 
point j,k lies just ahead or just outside of the leading edge, then the value of 
in equations (F5) and (F6) or the value of in equation (F7), respectively, must be 

the average of the values above and beneath the wing. If the point j,k lies at the trailing 
edge of the wing and the flow at this point is subsonic, the potential jump across the vortex 
sheet at the point j+l,k is set equal to that across the wing at the point j-l,k in order 
to satisfy the Kutta condition. 


The line relaxation process used here is the same as that used by Murman and Cole 
(ref. 8), Bailey (ref. 28), and others. A relaxation factor of 1 is used at all points where 
the flow is supersonic. At points where the flow is subsonic, the value of the relaxation 
factor is between 1 and 2. 
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SOLUTION PROCEDURE FOR GENERALIZED BROWN AND MICHAEL 
SEPARATED- FLOW MODEL FOR WING-BODY COMBINATIONS 


It appears that a computer program for the Brown and Michael model for separated 
flow past a general nonconical wing -body configuration has not been set up previously. 
Smith (ref. 49) has set up a computer program for nonconical wings, and Wei, Levinsky, 
and Su (ref. 31) have set up a program which uses the more complicated model of Mangier 
and Smith to calculate the flow over nonconical wing bodies. However, the wing bodies 
which can be treated by the method of reference 31 must have conical wing -body noses. 

The equations for the Mach number dependent Brown an^ Michael model for sepa- 
rated flow past a general nonconical wing -body configuration are 





|hyy 
2yy - vh 

y 

h 

y 

h (\2 * t /) 


v\ 




v^v ■ W 
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B. 


= (yy^ - 


2z 2 
V V 


B, 


= ^y,\{yv^ - 


E - 1 ’"'>*^y 


B 

B 2 + B ^ 
y z 


S " \®zO + ^z®yO 


~ ^y®yO " ^z^zo 


^ = iV - V ^0 - 2AyA^B,0 


^z 2Ay.A2.ByQ + ^Ay - A^ 

D , , - "v") 


°20 = 


6y z r, ^ 
V b 

(yy^ + 


„ _ ^yO^yO + ^zO^zO 

y *R 2 -u "R ^ 
®yO +®zO 
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For separated- leading- edge flow, 
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P -!v1£2 

y " Ri 


For a strake -generated vortex is constant and My, M^., Ny, N^, Oy, O^,, Py, 
and P_ vanish. 

2a 



Ry = Qy (yy\ + ^^v) - % (^v'^V - W) 
Rz " (^v^v “ ^v\) (^v\ ^v^v) 
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yy^yO + ^v^zO 


(yy^ + (®yO^ ^J) 


Gz = 


yy^zO - \^yO 


(yy^-^^^)(V-^®zO^) 


Hji = hRy + My 

Hi 2 = -hR^ + Ny 
Hji = hR^ + 

»22 = •'Ry + 


»1 = 


-sm a 


c F 

S 2 Uyh ^yj 



V + Z 
■^v V 


dr, 

2^^°y)‘'b-to 




H2 = sin a 


C + F 

z 2 Uyh 



y 2 + z 2 
Jy ^ y 


+ 2G„ r 


dr. 


^ dx 


- IVv^«zV-^Oz"''Gz)s-Pz^ 


dx 


H1H22 - H2H12 
^11^22 “ ^21*^12 
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dr 


dx 


H2H1I - H1H2I 

% 1^22 ■ ^ 21^12 


These equations are solved by fourth-order Runge-Kutta integration. The function Rj^ 
has a value of 1 if the attached-flow solution is obtained from slender-wing theory. If 
the approximate subsonic or supersonic attached-flow solutions are used, the function Rj 
is evaluated with equation (D41). Note that it is only through equation (D41) that the Mach 
number influences the solution. 

In order to start the solution at a point just downstream of where the leading edge 
of the wing emerges from the body surface, the configiiration in the vicinity of the junc- 
tion is approximated by a slender wing with the semispan y 2 - « r^^ at an angle of 

attack of 2o!. This approximation is depicted in figure 40. 


Exact Approximate 




Figure 40.- Approximation used to start calculation. 

The empirical formulation of Jobe (ref. 50) is used to determine initial values for 
y^ and z^. The starting equations are 

( R-i sin o\ /R-i sin 

-i 1 - 0.12152 

v-vj VVi’ 

Vv = Tb + (yj - i-b) - 0.98407zp + 2.29981Z|,^ - 1.78322Z|,®j 


= (y2 - 


7 = 



- --b)^ - (^2 - ■■b)^ 
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_ ^v(yy - ^b) 

r, 2 





- (y2 - + 4772^2^^^ 


= u 

\y 
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